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PREFACE 



Idcmpotcnt mathematics is a new branch of mathematical sciences, rapidly 
developing and gaining popularity over the last decade. It is closely related 
to mathematical physics. Tropical mathematics is a very important part of 
idempotent mathematics. The literature on the subject is vast and includes 
numerous books and an all but innumerable body of journal papers. 

An important stage of development of the subject was presented in the 
book Idempotency edited by J. Gunawardena (Publ. of the Newton Insti- 
tute, vol 11, Cambridge University Press, Cambridge, 1998). This book 
arose out of the well-known international workshop that was held in Bris- 
tol, England, in October 1994. 

The next stage of development of idempotent and tropical mathemat- 
ics was presented in the book Idempotent Mathematics and Mathematical 
Physics edited by G.L. Litvinov and VP. Maslov (Contemporary Math- 
ematics, vol. 377, American Mathematical Society, Providence, Rhode 
Island, 2005). The book arose out of the international workshop that was 
held in Vienna, Austria, in February 2003. 

The present volumes contain materials presented for the international 
workshop Idempotent and Tropical Mathematics and Problems of Math- 
ematical Physics (Moscow, Russia, August 25-30, 2007). 

It is our pleasure to thank the Independent University of Moscow and 
the Poncelet Laboratory of this university as well as the Russian Fund 
for Basic Research and CNRS (France) for their important support. We 
are grateful to a number of colleagues, especially to L. Kryukova and M. 
Tsfasman of the Poncelet Laboratory, T. Korobkova and Yu. Torkhov of 
the Independent University of Moscow, and A. Sobolevskh of the Moscow 
State University, for their great help. We thank all the authors of the 
volumes and members of our "idempotent /max-plus/ tropical community" 
for their contributions, help, and useful contacts. 

The editors 
Moscow, August 2007 
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Representation of stationary solutions of 
Hamilton-Jacobi-Bellman equations: a 
max-plus point of vie^\Q 

Marianne Akian 

This talk gathers two recent works: the first one is a joint work with 
S. Gaubert and C. Walsh first presented in [5], the second one is a joint 
work with B. David and S. Gaubert presented in pQ. 

1. Nonlinear eigenfunctions and stationary solutions of Hamil- 
ton-Jacobi-Bellman equations 

Let us consider a diffusion control model on a subset X of M. n , that is a 
stochastic process x t with values in X satisfying the stochastic differential 
equation 

(1.1) tfx t = g(x t , u t ) dt + cr(x t ,u t ) db t 

where (b t ) is a p-dimensional brownian motion, u := (u t )t>o (the control) 
is a stochastic process with values in a subset U of M. p and adapted to the 
filtration of (h t ), and the drift g : X x U — > K ra and the standard deviation 
a : X x U — ► M n , p (K) are given. 

The stochastic control problem with horizon T consists in maximizing 
over all the controls u the quantity 

i-T 



(1.2) E 



/ L(x s ,u s ) ds + 0(x T ) 
Jo 



where x t is the solution of (1.1 1 with initial condition xo = x, L is the 
Lagrangian and cf> is a final reward. Let us denote by v T (x) the value of 
this optimization problem, and by S T the map which associates v T to <fi. 
The familly of operators {S t }t>o is the (non linear) evolution semigroup 
associated to the control problem. Moreover, each operator S l is additively 
homogeneous (£*(/U + 0) = fi + S" (</>), where (u + 4>){x) = ^ + <j>(x)), and 
order preserving, thus it is nonexpansive for the sup-norm. If the control 
problem is purely deterministic then the operators S are max-plus linear, 
that is S t {4> V tp) — S 1 ^) V S l (ip). In general, the operators S* are convex, 
which means that, for all t > and x G X, the map i— > S t (c/))(x) is 



convex on 



We say that A is an additive eigenvalue of the evolution semigroup if 
there exists a function d> : X — + K such that for all i > 0, S 1 ^ = Xt+6. The 



: The present work was partially supported by the RFBR/CNRS grant 05-01-02807. 



Stationary solutions of HJB equations 



7 



function <f> is called an additive cigcnfunction of the evolution semigroup 
associated to A. If X is compact and we restrict ourselves to continuous 
eigenfunctions, the semigroup has at most one eigenvalue. Moreover, under 
some regularity assumptions on L,g,a, the eigenfunctions are exactly the 
viscosity solutions <fi of the crgodic Hamilton-Jacobi-Bellman equation 

(1.3) \-H(x,D(j>(x),D 2 <j>(x))=Q, xeX, 
where the Hamiltonian of the problem is given by 

(1.4) H(x,p, A) = max(- tr(cr(x, u)a(x, u) T A) + (p, g(x, u)) + L(x, u)). 

In that case, A is the maximal mean reward by unit of time (the ergodic 
reward) . 

Given an eigenvalue A of (<S')t>o, we are interested in characterizing 
the set S\ of associated eigenfunctions, or of solutions (j> of (1.3 1. 



2. Related results 

In the particular deterministic case (a = 0), the discrete-time analogue of 
this problem consists in the characterization of eigenvectors of max-plus 
linear operators, which has received a considerable amount of attention, 
see for instance |8j 1161 HI [3]. In the finite dimensional setting, eigen- 
vectors are max-plus linear combinations of extremal generators which are 
themselves in bijection with the "critical classes" (critical classes can be 
seen as the max-plus analogue of recurrent classes) . The deterministic con- 
tinuous time problem itself has been studied by Maslov, Kolokoltsov, Sam- 
borskii and other members of the "idempotent analysis" school |18l 119] , 
and by Rouy and Tourin |20j in some special cases. More recently, it has 
been studied as a part of the "weak KAM" theory developped by Fathi 
|14L I13L I12j and Fathi and Siconolfi |15j . In this setting, it is shown that 
when X is a Ricmannian manifold and the Lagrangian has smoothness 
and strict-convexity properties, an cigcnfunction is uniquely determined 
by its restriction to the "projected Aubry set" . This set can be thought of 
as a continuous analogue of the set of "critical states" of the finite dimen- 
sional max-plus spectral theory. The case of a non-compact state space X 
has been studied in different settings by Contreras and by Ishii and 
Mitake [i~7j . 

In Section [3] we present briefly the results of [5], where we showed that 
general representation results hold in the continuous-time setting, without 
any regularity assumption on the Lagrangian. These results were inspired 
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by the discrete-time theory developped in [3] and rely on a compactifica- 
tion of the state space X, which is the max-plus analogue of the Martin 
compactification in potential theory, and is similar to the compactification 
of metric spaces by horofunctions (generalised Busemann functions) . 

It is natural to ask whether an analogue of the weak KAM theory can 
be developped for stochastic control problems. Such an analogue does exist 
in the simpler finite state space and discrete time case. Indeed, it is shown 
in [2] that the additive eigenvectors are determined by their restriction to 
a subset of "critical states" obtained by taking exactly one element in each 
"critical class" . Here critical states are defined in terms of subdifferentials, 
and can be interpreted as follows: a state is critical if there is an optimal 
stationary randomised strategy for which it is recurrent, and two critical 
states are in the same critical class if they are in the same recurrence class 
for an optimal stationary randomised strategy. In the continuous time 
stochastic case, characterization results were obtained in the uniformly 
elliptic case and under various settings, by Bensoussan [10 , Akian, Sulem 
and Taksar [;6 , and Barles and Da Lio [9]: the eigenfunction is then unique 
up to an additive constant. 

In Section [4] we present briefly the results of pQ giving a description of 
the additive eigenspace similar to the one of [2], in the simplest degenerate 
case in which there is only a finite number of " singular points" x\ , . . . , Xk 
playing the role of critical states and classes. In the deterministic case 
similar results were obtained in |20l 118] . 



3. Hamilton-Jacobi equations on non-compact spaces 

In [5], we consider general time continuous semigroups (S t )t>o of max- 

plus linear operators with kernel. We assume that : R — > R can be 
written as S t g(x) = sup^g^ (S^ir,?/) + g{y))- for some function (x,y) i— > 
S t (x, y) £ R. This includes the case of the evolution semigroup associated 



to the deterministic optimal control problem with dynamics (1.1) with 



(7 = and criteria ( 1.2 1. Without loss of generality, we assume that A = 0, 
in which case an eigenfunction is called a harmonic function. Harmonic 
functions may take the — oo value, so that the set £q of harmonic functions 
is a semimodule over the max-plus semiring R max (recall that this is the 
set R U {— oo} endowed with max as addition and + as multiplication). 
We need the following assumptions : 
(Al) S*(x, y) :— sup t>0 ^'(x, y) is finite for all x and y in X. 
(A2) For all t > and x, y £ X, S"(x, y) = sup 7 {/(7)}, where the 
suprcmum is taken over all paths 7 : [0, t] — > X from x to y, and 
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where the reward I("f) is defined as 

n-l 

J( 7 ) := inf £ S^-^&O^+l)), 

with the infimum taken over all finite increasing sequences (U), 
i £ {0, . . . , n}, in [0, t] with to — and t n = t. 
The (max-plus) Martin kernel of the semigroup (S )t>a with respect 
to the basepoint b e X is defined by: 

K(x,y) = S'(x,y)-S m (b,y) . 

The (max-plus) Martin space ,.4t of (S 1 )t>o is the closure in the topology 
of pointwise convergence of the set Jif := {K(-,y) | y E X} C R . Any 
element of ^ is super-harmonic, which means that it satisfies S <fi < <f>. 
For all functions £ : X — > M max and for all 77 e we set: 

Vdv) : = limsup (S*(b,x) + £(x)) , 

and if £ € ^# we set: H(rj,£) := ^(n). The kernel -ff extends the kernel 
5* from X x X to ^# x up to a normalization, since 

H(K(-,x),K(;y)) = S*(b, x) + S*(x, y) - S*(b, y) . 

The minimal boundary of (S t )t>o is the set of elements £ of ^# that 



are harmonic and satisfy = 0. 



Theorem 3.1 ([5J Theorem 3.11]). f/nrfer Assumptions \(A1) - (A2)\ a 



function h : X — > f max is harmonic if and only if it can be written 
(3.1) h = sup (i'(iu) + w) , 

where v is some upper semicontinuous function from to R max - More- 
over, [ih is the greatest v satisfying this equation. 

Here, [ih plays the role of the spectral measure in the Martin repre- 
sentation theorem. In [5], we also show that the elements of the minimal 
Martin boundary are precisely the extremal generators of £q normalized 
in such a way that w(b) = 0, and that they are in correspondance with 
almost-geodesics . 

4. A degenerate Hamilton-Jacobi-Bellman equation on the torus 

In [1] , we study a simple degenerate stochastic control model on the torus 



T" = R"/Z n . We consider the stochastic model given by (1.1 1 with the 



criteria (1.2 1 on the set X = T n . We need the following assumptions: 
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(Al) L is C 2 on TP x U. 

(A2) g and a are Lipschitz continuous on T™ x U. 
(A3) L takes non-positive values. 

(A4) There exists k distinct points x\, . . . , Xk of T™ such that: 

(a) Vi = 1, . . . , k, 3ui G U such that g(xt, Ui) = 0, L(xi, Uj) = 
and a(xi,Ui) = 0. 

(b) Vx G T™ \ {xi, . . . , Xk}, Vu e J7, at least one of the following 
properties holds: (i) L(x,u) < or (ii) er(x, u)cr(x, u) T is a 
positive definite matrix. 

(A5) For all i — 1, . . . , k, there exists : T" — > [/ Lipschitz contin- 
uous satisfying u^(xi) = Ui, and a continuous function : 
T" — > R which vanishes on Xj and is positive elsewhere, and 
which is a viscosity solution of 

(41) - 1 -t r (a(x,u^(x)) < j(x,u^(x)) T D 2 W^(x)) 

-(.g(x,u W (x)),V^ w (x)) +L(x,u W (x)) > 0, x S T". 

Assumptions (A3) (A5) ensure that the point Xi is stabilized in probability 
by the control u^K Since X is compact, the evolution semigroup has a 
unique eigenvalue, and the previous assumptions imply that this eigenvalue 
is 0. The following result shows in particular that the set {xi , . . . , x^} plays 
a role analogous to the projected Aubry set. 

Theorem 4.1 QlJ). Under Assumptions (A1)\\(A5)\ the map £ Q — > 

WL k ,v i — ► (v(xi), ...,v(xk)) is a sup-norm isometry, whose image is a non- 
empty closed convex subset C ofM. k , that is invariant by all the translations 
(vi, ... ,Vk) !->■ (M + vi,... ,n + Vk) ([i. e R). 

A more precise description of the convex set C is given in pQ . 
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On the assignment problem for a countable state spaca^ 

M. Akian, S. Gaubert, and V.N. Kolokoltsov 

1. Introduction and formulation of the results. 

Our main results are formulated below as Theorems 1.1 - 1.3. The first two 
theorems are proved in Section 2. A proof of Theorem 1.3 will be given 
elsewhere. Section 3 contains the algebraic interpretation of our results 
and methods. 

Let X be either the set of natural numbers N or that of all integer 
numbers Z. Further we work with infinite matrices B — (bij), i,j € X that 
will always have entries in R U {— oo} and satisfy the following condition: 

(C) For any i there is a j such that bij =/= — oo, for any j there is an i 
such that b^ ^ — oo, and bij — ► — oo as \i — j\ — > oo. 

Any such matrix B defines the mapping B from the space of functions 
bounded below on X to the set M. x of all real valued functions on X, by 
the formula 

(1-1) (B/)i = BU P {6 ii -/ i } 

j 

By B T we shall denote the transpose matrix of B and the corresponding 
operator 

(1.2) (B T g)j = sup{fe iJ - <?i}. 

i 

A crucial fact about B T is that the pair (B, B T ) defines a Galois con- 
nection in M. x , which means in particular (see [2]) that B T is a generalized 
inverse to B in the sense that if the equation Bf = g with a given g £ ~R X 
has a solution / E R x , then necessarily / = B T g is also a solution of this 
equation. 

The infinite dimensional theory depends crucially on the class of func- 
tions, in which the solutions to the equation Bf = g are sought, and on 
the corresponding definitions of solutions to the assignment problem. Let 
loo] l\\ and denote respectively the spaces of functions s = (sj), i G X 
on X such that sup 4 |sj| < oo; sup n J2\i\< n l s «l < 00 5 an< ^ * ne (finite) limit 
limini^oo s n exists. Let I = l(X) be any of these spaces. 

Definition 1.1. A matrix B (satisfying (C)) will be called l-strongly 
regular if there exists a function g £ I such that (i) / = B T E I, (ii) / is 
the unique solution in M. x of the equation Bf — g and (iii) g is the unique 



iPartially supported by the joint RFBR/CNRS grant 05-01-02807. 
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solution in M. x of the equation / = B T g. In this case g (respectively, /) 
is said to belong to the l-simple image of B (respectively, of B T ). 

Of course, it follows from this definition that B is Z-strongly regular if 
and only if B T is strongly regular. 

Remark 1.1. One can show (though this is not obvious) that in the 
case of finite X our definition coincides with the standard definition of 
strong regularity given by P. Butkovic, sec [5J. In fact, we added a crucial 
additional condition in our definition, which turns out to be automatically 
fulfilled for finite, but not for infinite X. 

Definition 1.2. For any two bijections F, G : X i— > X, we define the 
distance between them by 

p(F, G) = sup \F(n) - G(n)\. 

n 

The binary relation F ~ G iff p(F, G) < oo is clearly an equivalence 
relation on the set of bijections defining the decomposition of this set into 
non-intersecting classes. We shall say that F is locally bounded if it is 
equivalent to the identity map. 

Definition 1.3. A bijection F : X i-s- X is called a (global) solution 
or a strong solution respectively to the assignment problem for a matrix 
B if 



(1.3) liminf \ (b iF{i) - b iG (i)) > 

\i\<n 



for any other bijection G : X i— > X or if 

(1.4) liminf V (b iF[i) - h lG[i) ) > 0, 

n — >oo * — ' 

\i\<n 

respectively. We say that this solution is a locally bounded l-solution, if F 
is locally bounded and the "solution sequence" bip ^) belongs to l(X). We 
say that F is a local solution if (1.3 1 (or (1.4 1 respectively) holds for all G 
such that the distance between F and G is finite. 

If a strong solution exists, then the solution to the assignment problem 
is obviously unique. 

Definition 1.4. A matrix B is called normal (respectively, strongly 
normal) if all its non-diagonal entries arc non-positive (respectively, nega- 
tive) and ba = for all i. 
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This definition is literally the same as the usual finite- dimensional 
one (see [4])- The normal (respectively, strongly normal) matrices present 
a class of examples, where the identity map is an obvious locally bounded 
Zi-solution (respectively, unique solution) to the assignment problem. As 
our first result will show, this class of matrices presents natural "normal 
forms" for strongly regular matrices. 

Definition 1.5. Matrices B and C are called (locally bounded) l- 
similar if there exist two locally bounded bijections H : X i— » X , K : X i— > 
X and two vectors g and / from l(X) such that 

(1-5) Cij = bjni) K (j)- (f>i-ipj. 

This is also a standard definition in the case of finite X (see e.g. [4j 
and Section 4 below for an intuitive interpretation). The importance of 
this notion is basically due to the following result. 

Proposition 1.1. (i) Conditions (C) and l-strong regularity for ma- 
trices with entries in R U {— oo} are all invariant under I -similarity, (ii) 
The property to have an If solution (in particular locally bounded or 
strong) to the assignment problem for matrices with entries in RU{— oo} is 
invariant under l\- similarity. (Hi) The property to have a locally bounded 
local Iq- solution (in particular strong) to the assignment problem for ma- 
trices with entries inlU {— oo} is invariant under Iq- similarity. 

Proof, (i) The invariance of condition (C) is obvious. The invariance 
of Z-strong regularity follows from the observation that if C and B are 



related by (1.5 1 then the equation g = Cf is equivalent to the equation 
(g + 0)(H- 1 ) = B((f + i/j)(K- 1 )). 

(ii) Let a bijection F be a solution to the assignment problem of a 
matrix C. Notice that 

(1.6) 

^ (CiF(i) - CiG(i)) = {bn(i)KF{i) - b H (i)KG{i)) + X! WteOO _ ^F(i))- 

Clearly, the last sum on the r.h.s. tends to zero as n — > oo whenever 
if> G It. Hence, F is an Zi-solution (respectively, a strong Zi-solution) to the 
assignment problem for the matrix C if and only if the mapping KFH^ 1 is 
an Zi-solution (respectively, a strong Zi-solution) to the assignment problem 
for the matrix B. 

(iii) By the previous discussion, it suffices to show that the sequence 

K = ^2 ^G{€) ~ i>F(i)) 

\i\<n 
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tends to zero as n — > oo if and only if both F and G are not infinitely 
far from the identity map. To this end, observe that due to the last 
condition, there exists a natural number p such that for every n the sets 
{F(i) : \i\ < n} and {G(i) : \i\ < n} both contain the set {i : \i\ < n — p}. 
Hence 

K = i>G(h) + V>G(* 2 ) + - + 4>G( lp ) - ^F( n ) - - - i>G(j P ), 

where ik (respectively jk) are such that |G(ij;)| > n — p (respectively 
l^(jfc)| > n — p). We recall now that the function ipn has a finite limit as 
n —> oo, which immediately implies that b n tends to zero (the statement 
that we wanted to prove) . This was the crucial application of this assump- 
tion, which seems to be the weakest possible to provide a link between the 
solutions to the assignment problem for similar matrices. □ 

Theorem 1.1. A matrix B (satisfying (C)) is l-strongly regular if and 
only if it is l(X) similar to a strongly normal matrix. 

It is of course interesting to know what can be said about the assign- 
ment problem for a regular matrix itself (not just for some of its similar 
matrices). In the analysis of this question (as well as the inverse one), an 
important role is played by the following functions describing in some sense 
the size of the problem. Namely, let Fbea (possibly local) solution to the 
assignment problem of a matrix B. Let "the optimal distance" between 
the points i,j be defined as 

^ 7 ^ hj = sup[fe iF(il) - b nF(il) + b llF{i2) - b l2F(i2) + ... 

+ b in _ 1 F(i n ) - b in F(i n ) + bi n F(j) " bjF(j)], 

where sup is taken over all n — 1,2,... and all collections ii, i^, i n of 
points from X, and "the potential" and "the inverse potential" as the 
functions on X are given respectively by 

(1.8) & = sup bij, fa = sup&ji. 

j j 

The following simple properties of these functions are crucial: 

(1) the values of b, <f> and 4> do no change if one takes the sup only 
over families i\ 1 ...,i n with pairwise disjoint points (in fact, any 
cycle gives a non-positive contribution due the assumption that 
F is a solution to the assignment problem) ; 

(2) bu, 4>i and fa are nonnegative for all i (in fact, take n = 1 and 



i\ = i in ( 1.7 1) 
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(3) the functions <f> and —(f> satisfy the equation 

(1.9) fi = sup[&ii?y) - b jF ^) +fj], VieX, 

j 

or, equivalcntly, 

(1.10) f = B1>, 1>j = b F -i U)j -f F -i j . VieX, 

(4) the functions <f> and —(f> satisfy the equation 

(l.H) /<=SU P [fcy + / i ] 1 ViGX 

j 

Observe that if B is normal then bij < and <fii — (pi — for all i, j. 

Theorem 1.2. (i) If a matrix B (satisfying (C)) is lo-strongly regular, 
then it has a (necessarily unique) locally bounded local strong l^-solution 
to its assignment problem such that 

(1.12) limsup&y < 0, 

and that the potentials <j) and cj> are bounded functions, (ii) If B is l\- 
strongly regular, then this solution is also a global l\-solution. 



To prove the converse to Theorem |1.2[ we shall use the following ad- 
ditional technical assumption on a solution to the assignment problem: 

(B(l(X))) Either the potential 4> or the inverse potential 4> belong to 
l(X). 

Theorem 1.3. If I is either l or l\ and if the assignment problem for 
a matrix B has a (possibly local) locally bounded strong l-solution satisfying 
condition (B(l(X))), then B is strongly l-regular. 

We have to indicate an unpleasant small gap between the necessary 
condition and the sufficient condition: from strong ^-regularity it follows 
that the potential <p belongs to loo , but in Theorem |1.3| we assume that 



4> G lo (which implies ( 1.12 ) ) . However, this discrepancy vanishes when 
we consider classes of similar matrices, as the following direct corollary of 
Theorem |1.1| and |1.3| suggests. 

Corollary 1.1. Let 1{X) be either Iq or l\. Then a matrix B is 
strongly l{X)-regular if and only if it is l(X)-similar to a matrix having a 
strong solution to the assignment problem satisfying condition (B(l{X))). 
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2. Coverings and sub-differentials. Proofs of Theorems and 

For the analysis of the equation Bf = g (also in a more general setting 
of uncountable X), an important role belongs to the notion of (abstract) 
sub-differentials. 

Definition 2.1. For a matrix B and /, g e R x , the abstract sub- 
differentials (or B-sub-differentials) are defined as follows (see [2] and ref- 
erences therein) 

df(j) = {keX: (Bf) k = sup{b kl - /,} = fty - /,•}, 

I 

d T g(i) = {keX : (B T g) k = sup{6 ;fe - ,g;} = b ik - g. t }. 

I 

For a given / the sub-differential is a mapping from X to the set 
P(X) of subsets of X. For any such mapping G its inverse mapping 
G- l Ih P(X) is naturally defined as G _1 0') = {i ■ 3 e 

We shall start with the following well known basic property of sub- 
differentials that we prove here for completeness. 

Proposition 2.1. If g = BB r g e R x , then (d T g)- 1 = dB T g. 

Proof. 

(5 T ff) _1 0') = {i : (B T g)j = sup{^ - = 6^ - 

i 

which is the same as gi — 6,j — (B T g)j, or equivalently B(B T g)i — bij — 
(B T g)j, and which means that i e d(B T g)(j). □ 

Proposition 2.2. Suppose that functions g, B T g € E x are bounded 
from below. Then B T g is a solution to the equation Bf — g if and only if 
<9 T <?(*) 7^ f or a M i or r equivalently, if the family of the sets {d T g) -1 (j) , 
j € X, is a covering of X. 

Proof. This is a direct consequence of a more general Theorem 3.5 
from [2], where one only has to observe that the assumption that / = B T g 
is bounded from below ensures that the set {j : bij — fj > (3} is finite for 
any i € X and [3 € K, which is the crucial condition for the applicability 
of this theorem. □ 

Definition 2.2. Let G be a mapping from X to the set of its subsets 
P(X) and let the family of subsets {G(j)}j & x be a covering of X. An 
element j G X is called essential (with respect to this covering) if 

3i e X : i i U keX \ : G(k). 
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The covering is called minimal if all elements of X are essential. 

Proposition 2.3. Suppose that functions g, B T g e R x are bounded 
from below. Then B T g is the unique solution of equation Bf — g in M. x 
if and only if {d T g) -1 (j) , j £ X, is a minimal covering of X. 

Proof. This is again a consequence of a more general Theorem 4.7 
from [2]. □ 

The key point in proving Theorems |1.1| and |1.2| is contained in the 
following statement. 

Proposition 2.4. Suppose that functions g, B T g £ R x are bounded 
from below, and such that f = B T g is the unique solution of equation 
Bf = g, and g is the unique solution to the equation B T g = f . Then there 
exists a locally bounded bijection F : X i— » X such that 

(2.1) j = F{i) <=► df(j) = {*} <=► d T g{i) = {j} <=► j = 

In particular, 

Vfc^F(i) b iF (i)- f F {i)>hk- fk, 
Vk^i hF{i) - 9i> hF(i) ~ 9k- 

Remark 2.1. As one easily checks, the inverse statement holds as 
well: if a locally bounded bijection F and bounded below functions /, g 



satisfy (2.1 1, then / = B g is the unique solution of the equation Bf 



and g is the unique solution of the equation B g = f . 



PROOF of Prop. 12.41 Applying Proposition |2.3| to the equation 
B T g — f one concludes that for all i there exists j such that j E (9/) _1 (i), 
but j ^ (9/)_i(fc) for any k ^ i. In other words (df)(j) = {i}, which by 
Proposition |2.l| means that (d T g)^ 1 ^) = {«}■ Hence, defining the map- 
ping F : X i— > P(X) by the formula 



= {j : (df)(j) = «} = {j : (d T g)- 1 (j) = {z}}, 
one deduces that F(i) ^ for all i and F is injective in the sense that 



F(i) n F{k) — whenever i ^ k. Applying now Proposition 2.3 to the 
equation Bf — g, one finds that for all j there exists i such that (d T g)(i) — 
{j}. From this, one easily concludes that each set F(i) contains precisely 
one point and that F is surjective, which finally implies that F is a bijection 



X i — > X such that ((2J) holds. 

Let us show that F is bounded. In fact, since B satisfies (C) and / 
is bounded from below, it follows that for any C > there exists iV such 
that bij — fj < —C whenever \i — j\ > N. On the other hand, as g is 
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bounded from below, 6jF(i) — /f(i) = 9% is bounded from below, and hence 
\i — F(i)\ < N for large enough N and all i. □ 



sition 



Proof of Theorem 11.11 . Let F be a bijection constructed in Propo- 
From the equation &jF(i) — /f(i) = Si h follows that {&iF(i)} £ 



2.4 



whenever f,g £ l(X). Hence the matrix C with entries 

Cij = b iF (j) - fp(j) - {biF(i) - fF(i)) 

is strongly normal and is ^-similar to B. □ 



PROOF of Theorem 11.21 . The existence of required solution fol- 
lows from Proposition 1 1.1 1 T his s olution is actually given by the bijection 
F constructed in Proposition 



2.4 



The boundedness of <fi and (f> follows of 
course from ( |1.11| . To prove the latter, one observes that according to the 
second inequality of (2.2 1 



hj < 9i-gj, 

and the r.h.s. of this inequality tends to as i, j — > oo, since g £ Iq. □ 
3. Algebraic interpretation 

A natural algebraic language for the analysis of discrete event systems and 
optimal control is supplied by the so called idempotent algebra, in partic- 
ular the (max, +)-algebra (see e.g. [6], [8]). The (max, +)-algebra deals 
with the semiring MU{— oo} equipped with the binary operations © = max 
and ® = + and with finite-dimensional semimodules over this semiring. 
The main impetus to the development of this algebra (and further its 
infinite-dimensional generalizations, see [8], [9]) was a simple observation 
that the basic Bellman operator 

(3.1) (5/)<=sup(&tf + /,•), i£X 

j 

of the optimal control theory is linear in this structure, i.e. 

B(ot! ®/i®a 2 ®/ 2 = Qi® B(fi) ®a 2 ® B{f 2 ) 

for ai,a 2 £ MU{— oo} and fx, / 2 £ (MU{— oo}) x ■ (Note that previously we 
denoted by Bf the operator which would now be denoted by B(—f); this 
was more convenient for the study of the inversion of B.) In fact the oper- 



ators of type (3.1 1 are the most natural (max, +)-linear operators, though 
they do not exhaust all of them (see e.g. pQ, [7], [2] and references therein 
for classical and recent results on this "kernel type" representations). 

The introduction of main notions and objects studied in this article was 
motivated by the development of the (max, +)-algebra that supplies a clear 
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intuitive interpretation for them. For instance, the strong regularity turns 
to be a linear algebraic problem connected with the inversion of matrices 
(or more generally linear operators having kernel representation). Our 
notion of similarity is obtained by rewriting the classical algebraic notion 
of similarity of matrices in (max, +). Next, the solution to the assignment 
problem 



sup V" b iF a) = ® F (g> ieX b iF u) 

turns out to be the (max, +) analogue of the classical algebraic notion of 
matrix permanent. Namely, solving the assignment problem means finding 
the (max, +)- permanent of a matrix (in our case infinite dimensional). 
If this solution is strong, then one says that this matrix has a strong 
permanent. 

An important tool in algebra is given by the so called Kleene star that 
for a given matrix A is defined by 

A* — ffi°° A k 

(the powers A k are understood in the operations of a given algebra). In 
(max, +)-algebra the elements of A* clearly define the longest path on the 
graph associated with A (see details e.g. in [3] for finite and respectively 
infinite space X) 7 and its columns are natural candidates for the solution of 
the eigenvalue - eigenvector equation for A (equation of type (1.9 1). Hence 
the non-surprising appearance of A* in our setting (our functions fry and 
potentials 4>i represent appropriately normalized elements of A*). As was 
observed in |10j . the functions fry turn out to be useful also in the anal- 
ysis of the Monge-Kantorovich mass transfer problem, which is a natural 
analogue of the assignment problem for general measurable (uncountable) 
state space X. 
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Dequantization of coadjoint orbits: the case of 
exponential Lie groups^ 

Alt Baklouti 

It is well known that the unitary dual G of an exponential solvable 
Lie group G is homeomorphic to the space of coadjoint orbits. Given a 
coadjoint orbit O, the orbit method enables us to construct the associated 
unitary and irreducible representation. The inverse procedure is called 
Dequantization and it consists in going backwards from an irreducible 
unitary representation tt : G — ► U (H) to the coadjoint orbit O and its 
associated geometric objects. Here, U (H) denotes the group of unitary op- 
erators on some complex inner product space H. Towards dequantization, 
we consider the Poisson characteristic variety of some topological unitary 
modules over a deformed algebra appropriately associated with the repre- 
sentation in question. In the case of nilpotcnt Lie groups, we showed that 
the Poisson characteristic variety coincides with the associated coadjoint 
orbit. For exponential Lie groups, we conjecture that such a variety coin- 
cides with the Zariski closure of the orbit. In this work, we prove such a 
conjecture for some restrictive classes of exponential Lie groups. 
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Quantum Pontryagin principle and quantum 
Hamilton-Jacobi-Bellman equation: a 
max-plus point of view 

Viacheslav P. Belavkin 

We exploit the separation of the filtering and control aspects of quantum 
feedback control to consider the optimal control as a classical stochastic 
problem on the space of quantum states. We derive the corresponding 
Hamilton-Jacobi-Bellman equations using the elementary arguments of 
classical control theory and show that this is equivalent to a Hamilton- 
Pontryagin setup. We show that, for cost functionals that are linear in the 
state, the theory yields the traditional Bellman equations treated so far 
in quantum feedback. A controlled qubit with a feedback is considered as 
example. 

1. Introduction 

Quantum measurement, by its very nature, leads always to partial infor- 
mation about a system in the sense that some quantities always remain 
uncertain, and due to this the measurement typically alters the prior to 
a posterior state in process. The Belavkin nondcmolition principle [3, 5j 
states that this state reduction can be effectively treated within a non- 
demolition scheme [5], [6] when measuring the system over time. Hence 
we may apply a quantum filter for either discrete [1] or time-continuous 
[3] non-demolition state estimation, and then consider feedback control 
based on the results of this filtering. The general theory of continuous- 
time nondcmolition estimation developed in [6] , [8] , [9] derives for quantum 
posterior states a stochastic filtering evolution equation not only for dif- 
fusive but also for counting measurements, however we will consider here 
the special case of Belavkin quantum state filtering equation based on a 
diffusion model described by a single white noise innovation, see e.g. [7]. 
Once the filtered dynamics is known, the optimal feedback control of the 
system may then be formulated as a distinct problem. 

The separation of the classical world from the quantum world is, of 
course, the most notoriously troublesome task faced in modern physics. 
At the very heart of this issue is the very different meanings we attach to 
the word state. What we want to exploit is the fact that the separation of 
the control from the filtering problem gives us just the required separation 
of classical from quantum features. By the quantum state we mean the 
von Neumann density matrix which yields all the (stochastic) information 
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available about the system at the current time - this we also take to be the 
state in the sense used in control engineering. All the quantum features are 
contained in this state, and the filtering equation it satisfies may then to be 
understood as classical stochastic differential equation which just happens 
to have solutions that are von Neumann density matrix valued stochastic 
processes. The ensuing problem of determining optimal control may then 
be viewed as a classical problem, albeit on the unfamiliar state space of 
von Neumann density matrices rather than the Euclidean spaces to which 
we are usually accustomed. Once we get accustomed to this setting, the 
problem of dynamical programming, Bellman's optimality principle, etc., 
can be formulated in much the same spirit as before. 



2. Notations and Facts 

The Hilbert space for our fixed quantum system will be a complex, sepa- 
rable Hilbert space f) . We shall use the following spaces of operators: 

A = 23 (f)) - the Banach algebra of bounded operators on rj; 

A* = T (f)) - the predual space of trace-class operators on f); 

S = & (rj) - the positive, unital-trace operators (states) on f); 

T = To (f)) - the tangent space of zero-trace operators on f); 

Tq = T (f))* - the cotangent space (see below). 

The space A* equipped with the trace norm \\g\\ 1 = tr|p| is a complex 
Banach space, the dual of which is identified with the algebra A with 
usual operator norm. The natural duality between the spaces A* and A 
is indicated by 

(2.1) ( e ,A):=ix{gA}, 

for each g G A+, A G A. The positive elements, in the sense of positive 
definiteness g > 0, form a cone T + of the real subspace T C A* of all 
Hermitian elements g — g\ and the unit trace elements g G T + normalized 
as WqWi — 1 are called normal states. Thus S = T+nT%, where 7i = 
{t G T : tr t = 1}, and the extremal elements g £ S of the convex set S C 
T + correspond to pure quantum states. Every state g can be parametrized 
as Q (?) — Qa ~ Q by a tangent element q £ Tq with respect to a given state 
go G S. We may use the duality ( |2.1[ ) to introduce cotangent elements 
p G Tq . Knowledge of (q,p) for each q G 7o will only serve to determine 
p G A up to an additive constant (as the q's are trace-free): for this reason 
we should think of cotangents elements p as equivalence classes 

(2.2) p [X] = {A G A : A = X + XI, for some Aet}. 
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The symmetric tensor power Afy m — A ® sym A of the algebra A is 
the subalgebra of 23 (t)® 2 ) of all bounded operators on the Hilbcrt product 
space f)® 2 = F)®f), commuting with the unitary involutive operator S = 
of permutations 771 <g> 772 1— > 772 ® 771 for any r/j e f). 

A map £(£,•) from .4 = 25 (h) to itself is said to be a Lindblad gener- 
ator if it takes the form 

(2.3) SS(t,X) = i[H{t),X]+Y J ZRAX), 

a 

(2.4) £ R (X) = R^XR - l -R)RX - l -XR^R 

with H self adjoint, the R a e A (and the summations in (3) understood 
to be ultraweakly convergent [27] for an infinite set {R a }). The generator 
is Hamiltonian if it just takes form i [H (t) , ■]. The pre-adjoint £' = £* 
of a generator £ is defined on the pre-adjoint space A* through the rela- 
tion (£' (g) ,X) = (q,Z(X)). We note that Lindblad generators have the 
property £ (/) = corresponding to conservation of the identity operator 
/ G A or, cquivalently, tr {£' (g)} — for all g G A*. 

In quantum control theory it is necessary to consider time-dependent 
generators £(i), through an integrable time dependence of the controlled 
Hamiltonian H{t), and, more generally, due to a square-intcgrable time 
dependence of the coupling operators R a (t) . We shall always assume that 
these intcgrability conditions, ensuring existence and uniqueness of the 
solution g (t) to the quantum state Master equation 

(2.5) j t e{t) = z'{t,e{t)) = ${t, e {t)), 

for all for t > to with given initial condition g(to) = go € S, arc fulfilled. 

Let F = F [•] be a (nonlinear) functional F [q] on A* (or on S C 
.4*), then we say it admits a (Frechet) derivative if there exists an A- valued 
function V g F [•] on A* (7^-valued functional on T ) such that 

(2.6) limi{F[. + H-F[-]}-(r,V e F[.]), 

h— >o n 

for each t e 4* (for each t e To). In the same spirit, a Hessian V® 2 = 
V e <E> V e can be defined as a mapping from the functional on S to the 
4® 2 m -valued functionals, via 

Jm q J L{F[. + hT + h'r'] - F [• + hr] - F [• + fcV] + F [■]} 

(2.7) = (r®T',V e ® V B F[-]). 
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and we say that the functional is twice continuously differentiable whenever 
V® 2 F [•] exists and is continuous in the trace norm topology. 

Likewise, a functional / : X h-> / [X] on A is said to admit an A*- 
derivative if there exists an A*- valued function V xf [•] on A such that 

(2-8) lim \ {.f [■ + M] - / [■]} = (V x f [■] , A) 

for each A e 05 (f)). The derivative V xf [•] has zero trace, Vx/ [-4] € % 
for each A e A, if and only if the functional / [X — XI] does not depend 
on A, i.e. is essentially a function / (p) of the class p [X] £ T Q *. 

With the customary abuses of differential notation, we have for in- 
stance 

V,/ ((q, X)) = f ((q, X))X, Vxf «f?, X)) = f ((q, X)) g, 

for any differentiable function / of the scalar x = (g,X). Typically, we 
shall use V e more often, and denote it by just V. 

3. Quantum Optimal Control 

From now on we will assume that the Hamiltonian H and therefore ■& 
(and v) are functions of a controlled parameter u G hi depending on t 
such that the time dependence of the generator £ is of the form £ (u (t)). 
Moreover, we do not require at this stage the linearity of ■& (u, g) with 
respect to g, as well as the quadratic dependence a (g), which means that 
what follows below is also applicable to more general quantum stochastic 
kinetic equations 

dg. (t) =$(u (t) , g. (t)) dt + a (g. (t)) dw (t) 

of Vlassov and Boltzmann type, with only the positivity and trace preser- 
vation requirements tr {$ (u, g)} = = tr {a (g)}. A choice of the control 
function {u (r) :r£ [to>*]} 1S required before we can solve the filtering 
equation (Bclavkin equation) at the time t for a given initial state g a at 
time to- From what we have said above, this is required to be a W-valued 
function which we take to be continuous for the moment. 

The cost for a control function {u(r)} over any time-interval [t, T] is 
random and taken to have the integral form 

(3.1) ^l{u(r)};t,g} = J C (u (r) , q u (r)) dr + G (q u (T)) 

where {g, (r) : r e [t, T]} is the solution to the filtering equation with ini- 
tial condition g, (t) = g. We assume that the cost density C and the 
terminal cost, or bequest function, G will be continuously differentiable 
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in each of its arguments. In fact, due to the statistical interpretation of 
quantum states, we should consider only the linear dependence 

(3.2) C(u,g) = (g,C(u)), G(g) = (g,G) 

of C and G on the state g as it was already suggested in |4] . [6] . [10] . We 
will explicitly consider this case later, but for the moment we will not use 
the linearity of C and G. We refer to C (u) £ A as cost observable for 
u £ hi and G £ A as the bequest observable. 

The feedback control u (t) is to be considered a random variable (t) 
adapted with respect to the innovation process w(t), in line with our 
causality requirement, and so we therefore consider the problem of min- 
imizing its average cost value with respect to {u, (t)}. To this end, we 
define the optimal average cost on the interval [t, T] to be 

(3.3) S(t,g) := inf E [J. [{«. (r)} ;t,g]] , 

{«.(r)j 

where the minimum is considered over all measurable adapted control 
strategies {it. (r) : r > t}. The aim of feedback control theory is then to 
find an optimal control strategy {u* (t)} and evaluate S (t, g) on a fixed 
time interval [io,T]- Obviously that the cost S(t,g) of the optimal feed- 
back control is in general smaller than the minimum of E [J. [{u} ; t, g]] 
over nonstochastic strategies {u(r)} only, which gives the solution of the 
open loop (without feedback) quantum control problem. In the case of 
the linear costs ( |3.2| this open-loop problem is equivalent to the follow- 
ing quantum deterministic optimization problem which can be tackled by 
the classical theory of optimal deterministic control in the corresponding 
Banach spaces. 



3.1. Bellman Si Hamilton-Pontryagin Optimality. Let us first con- 
sider nonstochastic quantum optimal control theory assuming that the 
state g(t) £ S obeys the master equation (2.51 where i) (u, g) is the ad- 
joint £' (u) of some Lindblad generator for each u with, say, the control 
being exercised in the Hamiltonian component i [•, H (u)] as before. (More 
generally, we could equally well consider a nonlinear quantum kinetic equa- 
tion.) The control strategy {u (t)} will be here non-random, as will be any 
specific cost J [{u} ; to, go}. As for S (t, g) — inf J [{u} ; t, g] at the times 
t < t + e < T, one has 

S (t, g) = 

inf J [ t+£ C (u (r) , g (r)) dr + f C (u (r) , g (r)) dr + G (g (T)) 1 . 



t+e 
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Suppose that {u* (r) : r G [t,T]} is an optimal control when starting 
in state g at time t, and denote by {g* (r) : r € [t, T]} the corresponding 
state trajectory starting at state g at time t. Bellman's optimality principle 
observes that the control {u* (r) : r G [t + e, T]} will then be optimal when 
starting from g* (t + e) at the later time t + e. It therefore follows that 



S (t, g) = inf 

Mr)} 



t+e 



C (u (r) , g (r)) dr + S{t + e,g(t + e))>. 



For e small we expect that g (t + e) = g + tf (u(t) , g) e + o(e) and provided 
that S is sufficiently smooth we may make the Taylor expansion 



(3.4) S(t + e,g(t + e)) = 



l + e^- t +e{d {u (i),e),V) 



S(t,g) + o(s). 



In addition, we approximate 

f-t+S 



C (u (r) , g (r)) dr = eC (u (t) , + o (e) 



and conclude that (note the convective derivative!) 



S (t, g) = inf 



S(t,g)^+o(e) 



where now the infimum is taken over the point-value of u (t) = u G U. In 
the limit e 0, one obtains the equation 



(3.5) 



-S (t, g) = inf {C (u, g) + (t? («, g) , VS (t, . 



where V = V e . The equation is then to be solved subject to the terminal 
condition 



(3.6) 



S(T,q) = G(q) 



We may introduce the Pontryagin Hamiltonian function on To x 7g* 
defined by the Legcndrc-Fcnchcl transform 

(3.7) H#(q,p[X}) := sup {(#(u,Q(q)), XI -X)-C(u,g(q))}. 

Here we use a paramctrization g (q) = go — q, q G Tq and the fact that the 
supremum does not depend on A G R since ($ (u, g) ,1) = 0. Therefore 
H depends on X only through the equivalence class p [X] G T * which is 
referred to as the co-state. It should be emphasized that these Hamilto- 
nians are purely classical devices which may be called supcr-Hamiltonians 
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to be distinguished from H. We may then rewrite (3.5 1 as the (backward) 
Hamilton- J acobi equation 

(3.8) - (i, g (q)) + Ho (q,p [VS (t, g)} (q)) = 0. 

Applying the derivative V g = —V to this equation to q = go — g in 
the tangent space 7q we obtain the dynamical equation p = — \7 q Ti.^ (q,p) 
for the co-state pt = Qt (T, s) of the operator-valued function X (t, g) — 
VS (f , g), where Qt(T,s) — p[X (t, g)} is the solution of this equation 
satisfying the terminal condition Qt (T, s) = s := p [G] with p [G {g)\ (q) — 
p[G (g (q))] for G = VG. We remark that, if u* (q,p(X)) is an optimal 
control maximizing 

K* (u, q,p(X)) - (0 («, g (q)) ,\I-X)-C(u,Q («?)) , 

then the corresponding state dynamical equation 4zQ = •& (u* (g, X) , g) in 
terms of its optimal solution q t = go — gt (to, go) corresponding to g ta = go 
can be written as q — V p 7i^ (q,p), noting that 

(3.9) VpWtf (q,p) = Vp/C,j (u* (q,p),q,p) = -ti (u* (q,p),g{q)) 

due to the stationarity condition J^K-o (u,q,p) = at u = u* . This 
forward equation with go = for g* (to) = go together with the co-state 
backward equation with px — p [G] = s is the canonical Hamiltonian 
system. Thus we may equivalently consider the Hamiltonian boundary 
value problem 



(3.10) 



q t - VpHv (q t ,Pt) = 0, q = 
pt + V q Htf (q t ,Pt) = 0, p T = s 



which we refer to as the Hamilton-Pontryagin problem, in direct analogy 
with the classical case. The solution to this problem defines the minimal 
cost as the path integral 

S(*o,£o)= / [(q r \Pr) -H(qr,Pr)]dr + G(Q(q T )). 

J to 

Thus the Pontryagin maximum principle for the quantum dynamical sys- 
tem is the observation that the optimal quantum control problem is equiv- 
alent to the Hamiltonian problem for state and co-state {q} and {p} re- 
spectively, leading to optimality K§ (u,q,p) < TL^ (q,p) with equality for 
u = u* (q,p) maximizing (u,q,p). 
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1. Introduction 

Totally positive matrices play an important role in different areas of mathe- 
matics, from differential equations to combinatorics. Studying parametriza- 
tions of canonical bases, Berenstein, Fomin and Zelevinsky [1 established 
the so-called Chamber ansatz for flag minors of matrices. This result 
relies on Pliicker relations between flag minors. In this talk, we study 
functions which satisfy tropical Pliicker relations. We consider two ap- 
proaches to tropicalization of Pliicker relations. One approach is based on 
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tropicalization of the so-called special Pliicker relation (3-term relation) 
after writing it as a subtraction-free expression. On this way we get the 
class of TP-functions (on Boolean cube 2^) which can be seen as the trop- 
icalization of flag minors. We show that such functions are determined 
by their restrictions to the interval family of subsets of N, and that the 
class of submodular TP-functions coincides with the class of submodular 
functions on the interval family. This might be seen as the tropicalization 
of the Chamber ansatz. Our proof is based on a construction of DMTP- 
functions via normal flows in weighted digraphs. A DMTP-function is a 
function that satisfies tropical Pliicker relations, and at this point we con- 
sider the tropicalization in the sense of the second approach, meaning that 
we tropicalize an algebraic formula with subtractions in the right hand 
side as an inequality. 

2. Flag minors and Pliicker relations 

For X an n x n matrix and I a subset of N = {1, . . . , n}, denote by Aj(X) 
the determinant of the submatrix located in the intersection of the first 
|/| rows and the columns indexed by I. These determinants are called 
flag minors of X. It is known (see Fulton and Harris [3], p. 235), that flag 
minors of a matrix X satisfy the Pliicker relations 

A^u/Aauj = y^(-l) d(J,J ''' J ^A y4u ( J \ i ) U .,A j 4 U (j\ J -) Ui , 

for any pairwise disjoint A, I, J c N, \I\ > \J\, any fixed j S J, and an 
appropriate function <i(J, J, 

The Pliicker relations with 2 > |7| > | J| > 1 are of special interest. 
They are given by 

A^tuifcA^uj = AAuij^Auk + A^ujfcA^uj 
with any A and i < j < k, and 

^AUik^Aujl = ^AUij^Aukl + A Aujk^AUil 

with any A and i < j < k < I. 

Definition 2.1. A function F : 2 N — * M is said to be P-function if F 
satisfies the above Pliicker relations for all A and i < j < k, that is 

(2.1) F(A U ik)F(A U j) = F(A U ij)F(A U k) + F(A U jk)F(A U i). 

Denote by VT the set of P-functions. 

We collected properties of P-functions in the following 
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Theorem 2.1. a) Any P- function satisfies all Pliicker relations. 

b) Let F : 2 N — > M^g be a P-function. Then there exists a unique 
upper-triangular N x N matrix X such that F(I) = A/(X). 

c) Let I C 2 N denote the interval family constituted from intervals 
{i, i + 1, . . . , j}, i < j, and let res/ : M 2 — * R 1 denote the re- 
striction map from 2 N to X. Then the mapping res/ is a bijection 
between the subspace PF (c R 2 ) and R x . 



3. Tropical Pliicker functions 

We consider two ways of tropicalization of P-functions. Firstly, we can 
tropicalize (2.1) as follows: 

(3.1) f{AUik) + f{AUj^=max(f{AUij)+f{AUk),f{AUjk)+f{AUi)), 
with any A and i < j < k. 

Definition 3.1. A function / : 2 N — > R is said to be a TP-function 
if F satisfies the tropicalization (3.1) of Pliicker relations for all A and 
i < j < k. 

Secondly, we can think of tropicalization in the form of inequality: 

(3.2) f(A U I) + f(A U J) < max(/(A U (/ \ i) U j) + f(A U ( J \ j) U %)), 

for any pairwise disjoint A, I, J, |/| > |J|, and any fixed j G J. 

Specializing this to tropicalization of the 3-term Pliicker relation given 
by (2.1), we have that, for all disjoint Ac N and {i,j, k} C N, 

the maximum is attained at least twice among the three values 
(3.3) 

a = f{AUik)+f{AUj),b = f{AUij)+f{AUk), c= f{AUjk) + f{AUi). 

Similarly, for a 4-term Pliicker relation, that is, for all disjoint A C N 
and {i,j, k, 1} C N, we have that 

the maximum is attained at least twice among the three values 
(3.4) 

x = f(AUik)+f(AUjl), y = f(AUij)+f(Aukl), z = f{AUjk)+f(AUil). 

Definition 3.2. A function / : 2 N -> E is called a DMTF '-function if 
/ satisfies (3.3 1 and (3.4 1. 

We will show in the following section that the class of TP-functions is 
a subclass of DMTP-functions. 

We have the following property of DMTP-functions, which is closely 
related to valuated matroids, see Dress and Wenzel [2]. 
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Theorem 3.1. A DMTP-function f satisfies all tropical Pliicker re- 



lations 1 3.2) 



4. Flows in digraphs and DMTP-functions 

Here we propose a construction of DMTP-functions. 

We deal with a digraph G = (V, E), a function c : E — ► K of weights 
on the edges, and disjoint subsets S, T C V. We assume that \S\ = \T\ 
and that T is ordered: T = (£i,t2> • • • > £|t|), and denote {ti, . . . ,tj} by Ti. 

For F C E and X C V, denote: 

the sets of edges in F leaving X and entering X by and by 

Sp(X), respectively; 

the set S^(X) U S F (X) by S F (X); 

the number |<5j(X)| - \6 F (X)\ by div F (X). 

Definition 4.1. Let us say that F c E is a normal flow from 5" C 5 

if 

f i v « e 5', 

div F (w) = < -1 Vtie T| S ,|, 

[ for the other v eV, 

where div(u) stands for div({w}). 

This gives rise to the following important function / = f c on 2 s : 
(4.1) f(S') := max{c(F), F is a normal flow from S'}, S' C 5. 



Theorem 4.1. / defined in (4-1) is a DMTP-function 



The claim that TP-functions constitute a subslass of DMTP-functions 
can be obtained as a consequence of the following 

Theorem 4.2. Let f : 2 N — * K &e a TP-function. Then there exists a 
planar digraph G = (V, E) and a weight function c : E — ► R., swc/i iftai / 



is defined by (4-.1) 



Remark 4.1. For a TP-function on the Boolean 2^, we can take the 
unique planar of the following from: V := 1 < z,j < n}, and two 

edges emanate from the vertex which terminate either in (i — l,j) 

or in (i,j + 1), respectively, if both terminate points are vertices. 

The following property is important, see Kamnitzer [4], for applica- 
tions to crystal bases construction via MV-polytopes: 

Theorem 4.3. A TP-function f : 2 N — > K zs submodular (that is 
f{A) + f(B) > f{AUB) + f(AnB), A, B ClN) if and only if the function 
r esi(/) «s submodular on I. 
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Degree one homogeneous minplus dynamic systems and 
traffic applications: Part ]Q 

N. Farhi, M. Goursat, and J. -P. Quadrat 

We show that car traffic on a town can be modeled using a Petri net 
extension where arcs have negative weights. The corresponding minplus 
dynamics is not linear but homogeneous of degree one. Possibly depending 
on the initial condition, homogeneous of degree 1 minplus systems may 
be periodic or have a chaotic behavior (to which corresponds a constant 
throughput) or may explode exponentially. In traffic systems, when this 
constant throughput exists it has the interpretation of the average car 
speed. In this first part we recall the derivation of the 1-homogeneous 
dynamics of traffic system and show the existence of such systems with 
chaotic behavior and a constant throughput. 

1. Introduction 

At macroscopical level, the traffic on a road has been studied from different 
points of view, for example : 

• The Lighthill-Whitham-Richards Model [6 is the more standard 
one 

d t p+d x q = , 

q = f(j>), 

where q(x,t) — denotes the flow at time t and position x on the 
road, p(x, t) = denotes density, / is a given function, called the 
fundamental traffic law. It plays for traffic the role of the perfect 
gas law for the fluid dynamics. 
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• The kinetic model (Prigogine-Herman [7]) gives the evolution of 
the density of particles p(t, x, v) as a function of t, x and v the 
speed of particle 

d t p + vd x p= C(p,p) , 
where C(p,p) is an interacting term in general quadratic in p. 

The second model is more costly in term of computation time and 
therefore not used in practice. The first one supposes the knowledge of 
the function /. This function usually comes from experimental studies, 
or from theoretical studies using simple microscopic model. Here, we will 
recall a way to derive a good approximation of this law from a simple 
minplus linear system based on a Petri net. 

The main purpose of this paper is to generalize this fundamental law 
to the 2D cases where roads have crossings. The original minplus linear 
model on a unique road cannot be generalized easily in term of Petri nets. 
We have proposed in a previous paper a way to solve the difficulty by 
using Petri nets with negative weights. The dynamics of these Petri nets 
can be written easily. Being not linear in minplus algebra any more, they 
are homogeneous of degree 1. We recall here the derivation of these 1- 
homogeneous dynamics. 

In the first part of this paper we show that we can compute the eigen- 
values for these 1-homogeneous system but that chaotic dynamics may 
appear. In the second part we discuss the phases appearing in the funda- 
mental diagram, obtained numerically, and describe new situations where 
we can prove that the system is periodic. 

2. Traffic on a circular road 

Let us recall the simplest model to derive the fundamental traffic law on 
a single road. The simplest way is to study the stationary regime on a 
circular road with a given number of vehicles and then to consider that 
this stationary regime is reached locally when the density is given on a 
standard road. We present two way to obtain this law : - by logical 
deduction from an exclusion process, - by computing the eigenvalue of a 
minplus system derived from a simple Petri net describing the road with 
the vehicles. 

2.1. Exclusion process modeling. Following [3] we can consider the 
dynamic system defined by the rule 10 — > 01 apply to a binary word 
describing the car positions on a road cut in section (each bit representing 
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Figure 1. A circular road. 



a section 1 meaning occupied and meaning free see II in Figure [TJ>. Let 
us take an example : 

mi = 1101001001, m 2 = 1010100101, m 3 =0101010011, 
to 4 = 1010101010, m 5 = 0101010101, 

Let us define : - the density p to be the number of vehicles n divided by 
number of places m : p = n/m, - the flow q(t) at time t to be the number of 
vehicles going one step forward at time t divided by the number of places. 
Then the fundamental traffic law gives the relation between q(t) and d. 

If P < 1/2 then, after a transient period, all the vehicle groups split 
off and then all the vehicles can move forward without other vehicles in 
the way, and we have : 



If P _ 1/2 then the free place groups split off after a finite time and move 
backward without other free place in the way. Then m — n vehicles move 
forward and we have 



q(t) = q = n/m = d . 



q = (m 



ri)/m = 1 — d . 



Therefore : 
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P 



1/2 1 



Figure 2. The fundamental traffic law. 

2.2. Event Graph modeling. Consider the Petri net given in III of 
Figure [l] which describes in a different way the same dynamics. In fact 
this Petri net is an event graph and therefore its dynamics is linear in 
minplus algebra. The number of vehicles entered in the place i before 
time k is denoted x\. The initial vehicle position is given by booleans ai 
with takes the value 1 when the cell contains a vehicle and otherwise. 
We use the notation a — 1 — a, then the dynamics is given by : 

x k i +1 = min{ai_i + x*_ x , a, + xt? +1 } , 

which can be written linearly in minplus algebra : 

x^ 1 — a^—xX^_i © a^x^i . 

This event graph has three kinds of elementary circuits : - the outside 
circuit with average mean n/m, - the inside circuit with average mean 
(to — n) /to, - the circuits corresponding to make some step forward and 
coming back, with average mean 1/2, Therefore its eigenvalue is 

q = min(n/m, (to — n)/m, 1/2) = min(p, 1 — p) , 

which gives the average speed as a function of the car density. 

3. 2D traffic 

Let us generalize the second approach to derive the fundamental diagram 
to a regular town describe in Figure [3] The complete town can be modeled 
as a set of subsystems corresponding to a unique crossing and two adjacent 
roads. To write the dynamics of the town we have first to give the Petri 
net describing a crossing. 

A first trial is to consider the Petri net given in Figure [4] This Petri 
net is not anymore an event graph but following L. Libeaut[5] it is possible 
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Figure 3. A town. 



to write the nonlinear implicit minplus equation describing a general Petri 
net. In the case where the multipliers are all equal to one it is : 



(3.1) min 



+ x'{k-l)- =0, Vx : Vfc. 

where x{k) denotes the firing number of transition x and p a place of the 
Petri Net. But these equations do not determine completely the dynamics 




Figure 4. A simplified crossing. 

since solution to the Cauchy problem is not unique, in general. Indeed : - 
at place a n we may have a routing policy giving the proportion of cars going 
towards 2/2 and the proportion going towards 2/3 (which is not described 
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by the Petri net |4J) — at place a„ we may follow the first arrived the first 
served rule with the right priority if two cars arrive simultaneously at the 
crossing (which is also not described by the Petri net [3} . 

Precising the dynamics of Petri net in such way that the trajectories 
are uniquely defined corresponds to give another Petri net having only one 
arc leaving each place. Let us discuss more precisely these points on a 
simple system given in the first picture of Figure [5] 




Figure 5. Dynamic Completion. 



The incomplete dynamics of this system can be written in minplus 
algebra X ^ Xq — CLX . Clearly x^ and xz are not defined uniquely. 

We can complete the dynamics, for example, in the two following ways 
useful for the traffic: application : by precising the routing policy 

X ^ — X^ — *\J CLX 

- by choosing a priority rule 

X^ — (XX ~y / ^4 

~n _ r l T n - 1 r r 'n.-l In 

In the two cases we obtain a degree one homogeneous minplus system. 

This method can be applied to the crossing and we obtain a Petri net 
with negative weights which has only one arc leaving each place (that we 
call deterministic Petri net) see Figure [6] 



Tropical Pliickcr functions 



39 




FIGURE 6. A Complete Crossing. 



Neglecting the roundings the system can be written with minplus no- 
tations : 

Xi/S = di-iXi-i © OjXj+i, 

x nl 8 — Q'n x l x n+l/ x 2n © a n— l x n—l > 
< X 2n /S — a 2n XiX n+ i/(x n /S) © a 2n -lX2n-l , 

xi/5 = a n ^/x n x 2n © aix 2 , 

where <5 denotes the forward shifting operator acting on sequences. It is a 
general degree 1 homogeneous minplus system. 

Simulation of this system starting from shows that 

Yaax^/k = A, Mi . 

The constant A has the interpretation of the average speed. The fun- 
damental diagram gives the relation between the average speed and the 
vehicle density of the system. In Figure [7] we give this law in the cases of 
two circular roads with one crossing for different relative size of the two 
roads. We see that three phases appear on each fundamental diagram. 
These phases will be discussed in the second part of this paper. The 
experimental existence of this A motivates the study of the eigenvalue of 
1-homogeneous minplus system. 
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FIGURE 7. 2D-traffic fundamental diagrams. 

4. Eigenvalues of 1-homogeneous minplus systems 

The eigenvalue problem for 1-homogeneous system / : R™ in i— > R^in can 
be formulated as finding x € IR^in non zero i an d A £ JR min such that : 

Ax = /(x) . 

Since / is 1-homogcncous, supposing without loss of generality that if x 
exists x\ ^ e, the eigenvalue problem becomes : 

'A = h{x/x{) , 

< X2/X1 = {h/h){x) , 

x n jx Y = (f„/fi)(x) , 

Denoting y = {x^/x-i, • • • , x n jx\) the problem is reduced to the computa- 
tion of the fixed point problem y = g(y) (with gi-\{y) = (/i//i)(0, y)) to 
compute a normalized eigenvector from which the eigenvalue is deduced 
by : A = /i(0, y). But now g is a general minplus function. 

The fixed point problem has not always a solution. There are cases 
where we are able to solve the problem - / is affinc in standard algebra, - 
/ is minplus linear, - / is positive power function. In the first case there 
is a unique eigenvalue as soon as dim(ker(/ — /)) = 1. 

In the two last cases, the problem can be reduced to the minimization 
of the average cost by time unit using dynamic programming methods. 
The corresponding fixed points are unique and stable. 

Moreover, since max(x,y) = xy/(x y) games problem are also 1- 
homogeneous minplus systems and the solution of the corresponding eigen- 
value problem is known. 
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In the general case we may have unstable fixed points that, never- 
theless, we can compute by Newton method (which is exactly the policy 
iteration) but which don't give the information about the asymptotic be- 
havior of the system anymore. In this case the asymptotic is obtained by 
an averaging based on invariant measure which may be difficult to com- 
pute. Let us give an example of chaotic system which has a 1-homogeneous 
minplus dynamics. 

5. A Chaotic system example 

Let us consider the 1-homogeneous minplus dynamic system 

^ +1 = (4) 2 /4®2(4) 3 /(^) 2 , 

The corresponding eigenvalue problem is 




Figure 8. Cycles of tent transformation. 

The solutions are A = and y = X\/x2 satisfying the equation 

y = y 2 ®2/y 2 1 

which has for solutions y = and y — 2/3. These two solutions are 
unstable fixed points of the transformation f(y) — y 2 2/y 2 . But the 
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system y n +\ = f(yn) is a chaotic system since / is the tent transform 
(see [2] for example for a clear discussion of this dynamics) . In Figure [8] 
we show the graph of /, / o / , / o f o /, their fixed points and the 
corresponding periodic trajectories. 

In Figure [9] we show a trajectory for an initial condition chosen ran- 
domly with the uniform law on the set {(i — 1)/10 5 , i — 1, • • • , 10 5 }. The 
diagonal line in the picture is a decreasing sort applied to the trajectory. It 
shows that the invariant empirical density is uniform. We can prove that 




Figure 9. A tent iteration trajectory. 

the tent iteration has a unique invariant measure absolutely continuous 
with respect to the Lebesgue measure : the uniform law on [0,1]. 

More generally a chaotic 1-homogeneous minplus system will grow 
linearly with a value A given by : 




h(y)dKv) , 



where /i is the invariant probability measure of y depending on the initial 
value y°. For example, according to the initial value y°, the tent iterations 
y k stay in circuits or follow trajectories without circuit (possibly dense in 
[0,1]). 
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Degree one homogeneous minplus dynamic systems and 
traffic applications : Part 

N. Farhi, M. Goursat, and J. -P. Quadrat 

In this second part we discuss the phases appearing in the fundamental 
diagrams of traffic systems modeled by 1-homogeneous minplus dynamics 
and show the improvement obtained by traffic light control. 

We have shown in the first part that 1-homogeneous systems may 
have a chaotic behavior. Here we give a new subclass of 1-homogeneous 
dynamics having periodic trajectories. It generalizes the standard cases 
which need a monotony property. Moreover we show that this new, but 
still restrictive class, has applications to regular town traffic with crossings 
but without turning possibilities. 

1. The traffic fundamental diagram phases. 

The fundamental diagrams of quite different systems are similar to the 
one given in part I. We have studied the cases of two circular roads with 
one crossing and two crossings and the cases of regular towns with various 
number of roads on a torus. In all these cases we suppose the existence of 
right priority. 

The fundamental diagrams have always three phases corresponding 
respectively to low, average and high densities. We see on the fundamental 
diagram of Part I that : - for low densities the flow increases linearly with 
the density, - for average densities the flow is constant, - for high density 
there are deadlocks and the flows are null. 

On Figures (JlJ , ^ and ^ we show the typic asymptotic distribution 
of vehicles in the three phases [3j |4] . 
We see that : 
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Figure 1 . Two circular roads with one crossing case. Car 
distributions in the low average and high density phases. 




Figure 2. Four roads with two crossings. Car distribu- 
tions in the low average and high density phases. 




Figure 3. A regular town. Car distributions in the low 
average and high density phases. 

• Low density phase. There are so few vehicles in the network 
that after a transient regime, they move without obstructing each 
other on the roads and in the crossings. Thus, the "priority to 
the right" is not used, the vehicles moves as on a unique circular 
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road and the average flow is equal to the vehicle density in the 
network. This phase corresponds to densities less than 1 /4. 

• Average density phase. When the density is between 1/4 and 
1/2 (in the symmetric road cases), the vehicles can neither move 
freely on the roads, nor avoid each other on the crossings. There- 
fore "priority to the right" happens. The car on the priority road 
move freely and the waiting cars are all in the non priority road. 
The flow reaches the maximum value 1/4 corresponding to the 
full use of the crossings. 

• High density phase. When the density exceeds a quantity equal 
to 1/2 (in the symmetric case), at asymptotic regime, a closed 
circuit of vehicles on some nonpriority roads appears which cre- 
ates a complete deadlock of the system. 



2. Traffic light control. 

To avoid the deadlock due to right priority we can use traffic light controls. 
A Petri net describing the junction with the traffic light control is shown 
on Figure [4] The negative weight extension of Petri net is necessary to 
model the light phases in a time invariant way. The part of the Petri 
modeling the light control corresponds to the places a gi a c ,a gi a c . As long 
as a c contains a token a c = 1 the green light is for the North street, when 
a c = 1 the green light is for the East street. As long as a c = 1 we have 
a g = 1 and q v is authorized to fire (since thanks to the loop q v ,a g ,q v as 
soon as a token is consumed another one is generated in the place a g ). 
The main point is that when the token in a c goes in a c (phase change) the 
tokens in a g must be removed (this is done by the input arc with weight -1 
of the place a g ). More generally without negative weight we cannot model 
tokens staying less then a prescribed time. 

In Figure [5] we compare the fundamental diagrams of three crossing 
policies for a system composed of two circular roads of same size with two 
junctions. The three policies are : - right priority, - standard given phase 
duration, - feedback controlled duration (based on the road congestion) 
computed by LQG method. 

The control improves the average and the high density phases, without 
spoiling the low density one. The improvement given by the feedback 
control achieves the throughput obtained on a unique circular road without 
crossing but doubling the time spent in a place representing the crossing 
place. 
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Figure 4. Traffic lights modeling. 

Furthermore, the feedback control dissolves more efficiently the jams 
(that can appear locally in transient regimes) than the other policies would 
do. 




Figure 5. Comparison of three policies of managing the 
crossings : -right priority to the right (1), -open loop light 
control (2), -feedback light control (3). 



3. A subclass of triangular homogeneous dynamics 

In this section we study a subclass of 1-homogeneous minplus linear sys- 
tems for which we can prove the periodicity. Their dynamics belongs to a 
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where {u/c}fceN and {^fcjfceN are minplus column vectors, C is a minplus 
square matrix, A(uk) and B{uk) are two minplus O-homegeneous matrices 
depending of Uk- We call this type of systems Triangular 1-Homogeneous 



We call linear periodic dynamic (LP) a dynamic given by : 

x k +i = A k ®x kl xq given, 

where Ak are minplus matrices periodic in k. 

We can prove the following theorems (see the proofs in [5]). 

Theorem 3.1. Every T1H dynamics behaves asymptotically as a LP 
dynamics. 

Theorem 3.2. A T1H system with A{u) irreducible for every u € 
^■min satisfies : 



where : - Uq denotes the initial condition of u, - [i x (u ) = lim^^ Xk/k, 
- V is the set of the minplus eigen vectors of the matrix C . 

Theorem 3.3. Every LP dynamic yk+i = Ek <8> yk, such that the 
matrices Ek have the same support, is realizable by a T1H dynamics. 

4. Application to traffic 

We show that the traffic of regular towns with traffic light, buffered junc- 
tion but without turning possibilities can be modeled with a T1H dynam- 
ics. 

On the Petri net of Figure [6] the traffic light is modeled by the sub- 
system corresponding to the transitions u%, U2, U3, U4, which has no input 
coming from the rest of the system. The dynamic of this subsystem is min- 
plus linear. If the initial condition uq — (0, 0, 0, 0) the number of tokens in 
the places ao and b is boolean and periodic. To a cycle corresponds the 
four phases given in the Table [l] The junction has a buffer place in each 
direction (ai,6i) to avoid blocking. The phases 2 and 4 gives the time, 
for car entering in the junction, to go in the buffer and then to free the 
crossing. Indeed, a vehicle entering in the crossing (represented by the two 
places a n and b m ) leaves it surely in one unit of time. 



(T1H). 



max fix(uo) 
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Figure 6. Traffic light intersection without possibility of turning. 



Phase 


a 


bo 


Vertical light color 


Horizontal light color 


1 


1 





green 


red 


2 








red 


red 


3 





1 


red 


green 


4 








red 


red 



Table 1. The phases of the traffic light. 



The green duration of phase 1 and 3 is the sojourn time of tokens in 
the place ifi. The phases 2 and 4 have a duration of one unit. 

Proposition 4.1. The dynamics of the Petri net of Figure [5] has the 
T1H dynamics : 



,k+l 



'r2 



Y?3 



- x k+r 




'A x {u k ) 




x k 


z k+i 




A 2 {u k )_ 




z k 
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where ■ denotes oo, with 

ja u 1 /u 2 if = (n,n), 



and 



| independent ofu elsewhere. 

\b a u 3 /u 4 if(i,j) = (m,m), 
) independent of u elsewhere. 



We are able to explicit the asymptotic flows which are different ac- 
cording the direction followed by the vehicles. 

Theorem 4.1. The average flow on the horizontal (resp. vertical) 
road is given by A/4 where A is the unique eigenvalue of the irreducible 
matrix ^) 3 k=0 Ax(u k ) [resp. ® 3 k=Q A 2 (u k )]. 
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Max-plus cones and semimodulesj 

F. Faye, M. Thiam, L. Truffet, and E. Wagneur 

The concept of moduloi'd over a dioi'd has been introduced in \7[ . These 
algebraic structures have been considered in the context of production 
systems [5] , computer systems [9] , network systems [4] , or more generally 
for the modeling and analysis of discrete event systems pQ, [§]. For G. 
Gondran and M. Minoux (0)> a moduloi'd over a dioid is the algebraic 
structure left invariant under the action of a matrix A with entries in a 
dioi'd (the "space" of proper "vectors" of A) . This structure is also very 
similar to that of band-space over a belt of R.A. Cunninghame-Green [6]. 

The problem of solving linear equations of the type Ax = Bx in the 
max-algebra has been considered by many authors (cf [4], and [2], where 
additional references may be found). In [2] the authors show how to com- 
pute all solutions to a system of linear equations over a totally ordered 
idcmpotent semificld. 

In [3] the authors consider subsets of the positive cone E™ endowed 
with the max operator as the first composition law, together with "+" as 
the second composition law. They show how to relate subsets of K™ to 
the concept of generating vectors and bases defined in [llj . 

Since the early years, the terminology evolved, and the concept of 
idempotent semimodule over an idempotent semiring (or semifield) has 
emerged as the counterpart of that of vector space or, more generally, of 
module over a ring. 

The most general definition of a (finite dimensional) idempotent semi- 
module M is the following one: take two matrices A t B of finite size with 
coefficients in a semificld S, consider the inequalities Ax < Bx, and then 
define M as the set of solutions to this set of inequalities. Another way to 
define (or to represent) a (finite dimensional) semimodule over S is to give 
its basis, for example as the (independent) columns of a matrix with coef- 
ficients in S, It is then natural to ask how to get from one representation 
to the other. 

The aim of this paper is to study n-dimensional semimodules over a 
completely ordered and complete idempotent semifield defined by a pair 
(AB). 

Recall that an idempotent semigroup (S, V) is ordered by the relation 
s, t € S , s < t <=> s V t = t. If (S, V,0) is an idempotent semigroup 
with the neutral clement 0, then is the least clement of S, since for 
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every s € S,sVO = s. We will assume here that the set of scalars S 
is an idempotent completely ordered semifield, which is complete, i.e. S 
is totally ordered and complete, with least clement 0, endowed with two 
composition laws: V, and • such that : 

i) (S 1 , V) is an idempotent commutative monoid, with neutral clement 0. 

ii) (S\{ 0}, •) is an abelian group — hence (S \ {0}, •, <) is an £-group -, 
with neutral element written 1. 

iii) • is distributive over V, 

iv) is absorbent (for every s G S, ■ s = 0) . 

An idempotent semifield is also called a dioid pQ. 

Next we give a brief summary of our talk. In Part 1, we state some 
general results on semimodulcs defined by a pair {A,B). In particular, 
we give an explicit formula for the semimodulc of solutions to a single 
equation ciiX < biX in terms of the coefficients ay, and by of cij and 
hi (here and in the sequel, scalar product is denoted by concatenation). 
Take a permutation <jj from the symmetrical group S n such that, for j — 
1, . . . , k , a>i,<7i(j) < h,(nU)' while b ii<Ti y) < a ii<Ti (j) , j = k + 1, . . . n. Let 
J k(i) = < j < ri)\ a i)CTj (j) < b. h<Jz(j )}. Clearly for every j g J k (i), 

&<ji{j) is a solution to a^x < fyx, where the e^'s are the elements of the 
canonical basis of S n . 

Also, for every j £ Jfe(i) , I Jk(i), we have: 

=K*i{f) = b i,c t (j)( 1 V ^"iW^w) = b <,<ri(j) V 6 i.^(0( a »> 4 (<) 6 i.o-<0'))- 

Hence e CTj (j) V (^a-m^wO')) 6 ^^) a ^ so satisfies < ^i^- 
We have proved the following statement. 

Proposition 1. Ifbijt a$, i/ien i/ie sei of solutions to a t x < bix is 
a semimodule Mi generated by the k (n + 1 — k) vectors given, for every 
j E J k (j), by: 

e ai ( 3 ), and e a .(j) V (a7^ w b iiCT . (j -))e CT . w , I £ J k {j). 

The semimodule M, is generated by the columns of Vi given by the 
concatenation over J k {i) of the matrices 

e <7i(j) v ( a t> < (i 2 ) & <^<(i)) e o-i(^)l ' ' ' l e ^(j) V (a i > j (^_ fe )^,<r i y))e CTi (c_ fc) ]> 

where , 4 , . . . , £ n _ k } = {i\ 1$ J k (i) } . 

In Part 2, we give a geometric interpretation of the results of Part 1. 
Let M p stand for the semimodule generated by the solutions to a p x < b p x, 
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with <7p e S n , and J q (p) = {r (1 < r < n)|a CTp ( r ) < &cr p (r)}- We have the 
following statement. 

Theorem 1. For 1 < k < n - I, we have Mj n M p ^ {0} ijff one o/ 
i/ie following conditions holds 

»M(i)fK(p)^_0 

»J V S(<)CpW - A °p,i i y) 6 i.^Cj)- 



«eJ 9 (p) 



jeJk(i) 



In Part 3, we solve the system of inequalities < Bx using combi- 
natorial decomposition of the inequalities OjX < biX. 

In particular, we show that aix < biX is equivalent to a series of inequali- 
ties : 

o" 1 (a a x a V (a i2 V & i2 )z2 V ... V (a in V 6 in )a; n ) < x 1 or 
o^ 1 ((oji V ba)xi V a i2 a; 2 V (a i3 V b i3 )x 3 V ... V (a in V b in )x n ) < x 2 or 
or 

Kn{( a H V V ... V a in X n ) < X n . 

The first inequality may be written in matrix form (add the trivial 
inequalities x 2 < x 2 , . . . , x n < x n ) as: Pnx < x, with 

b^au b^ 1 (a i2 V b l2 ) ■ ■ o^ 1 (a in V b in ) 
10 
1 



Pa 












It is well-known that Pnx < x P*\X = x - This equation has a 

nontrivial solution iff an < bn, and in this case the solutions are given by 
the columns of P* x . 

Proceeding similarly for each line, we get P 2 * , . . . , P^ , and all solutions 
lie in the concatenation [Pii|Pj2| ■ • • |Pm]- 

Then we look at all the intersections of the type Pit H Pje, etc. in 
a combinatorial way, and devise an algorithm for the solution to these 
systems of inequalities. 

Finally, we give a complete description, both algebraic and geometric, 
for cases n = 2, 3. 
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Duality of cluster varieties 

V. V. Fock and A.B. Goncharov 

Cluster variety is an algebraic variety (strictly speaking, a scheme) 
defined by combinatorial data by explicit set of coordinate charts and 
transition functions. More precisely, for any collection of combinatorial 
data, called seed one associates three varieties _4|i| , X\j\ , and T)\i\ . These 
varieties possess canonical pre-symplectic, Poisson and symplectic struc- 
tures, respectively. One defines also a discrete group acting on all the 
three types of varieties and preserving the respective structures. The man- 
ifolds <¥|i| and T>\i\ admit a quantisation (noncommutative deformation of 
the algebra of functions) which is also 2) |i| -invariant. 

Varieties admitting cluster descriptions are simple Lie groups, mod- 
uli spaces of Stokes parameters, moduli of flat connections on Riemann 
surfaces, configuration spaces of flags, TcichmuUcr spaces and their gen- 
eralisations, the spaces of measured laminations and some others. One of 
the important features of cluster varieties is that they are defined not only 
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over a field but also over semifields (semigroups w.r.t. addition and groups 
w.r.t. the multiplication). For example, one can consider Tcichmullcr 
space, space of measured laminations and the space of flat PSL(2,¥)- 
connections over a surface S as the same cluster manifold but defined over 
the scmificld M>o of positive real numbers, tropical scmificld M* (which is 
ordinary R as a set with maximum for the addition operation and ordinary 
addition for the multiplication), and a field F, respectively. 
Let us give the precise definitions: 

A cluster seed, or just seed, I is a quadruple (I, Iq, e, d), where 

i) I is a finite set; 

ii) Jo C I is its subset; 

iii) e is a matrix e^, where i,j e I, such that <G Z unless i,j e I . 

iv) d = {di}, where i £ I, is a set of positive integers, such that the 
matrix Sij — £ijdj is skew-symmetric. 

The elements of the set I are called vertices, the elements of I are 
called frozen vertices. The matrix e is called exchange matrix, the numbers 
{di} are called multipliers, and the function d on I whose value at i is di 
is called multiplier function. We omit {di} if all of them are equal to one, 
and therefore the matrix e is skew-symmetric, and we omit the set Iq if it 
is empty. 

An isomorphism a between two seeds is a map I = (I,Io,e,d) and 
I' = (I',I ,e',d') is an isomorphism of finite sets a : I — ► I' such that 
cr(/ ) = I'o' ^cr(i) = di and e auaj = e^. Observe that the automorphism 
group of a seed may be nontrivial. 

For a seed I we associate a torus X\ = (F x ) / , called X -torus, another 
torus Ai = (F x ) / , called X -torus and the third one V\ = (F x ) /X/ called 
V-torus or a double torus. We denote the standard coordinates on these 
tori by {xi\i <G /}, {ai\i £ 1} and {yi,bi\i G /}, respectively. 

The A^torus is equipped with the Poisson structure 

(1) {x^^Cj} — CijXiXj 

The ,4-torus is equipped with the pre-symplectic structure (closed 2- 
form to possibly degenerate) 



(2) 




The D-torus is equipped with the symplectic form 
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The inverse of this form is a nondegenerate Poisson structure which can 
be written as 

(4) {VuVj} = eijViVj, {Vubj} = tfcPyibj, = 

Observe that these sructures are constant in logarithmic coordinates. 

Isomorphism between two X-tori Xj and X\* is a map given in coor- 
dinates by Xcr(i) = Xi, where a is an isomorphism of the seeds. Observe 
that there are much less isomorphisms of A"-tori then just isomorphisms of 
the corresponding Poisson manifolds. Isomorphisms of A- and 2?-tori are 
defined analogously. 

There exist the following maps between the tori: 

(5) Aj^X I} x, \[a : 



(6) Ai x Ai -> T>i, y l = Y]_ Oj' 1 , h = ai/ai, 

3 

Here aj are coordinates on the second A\-f actor. 

(7) Vi-tX-i, Xi = yu 
and 

(8) Vi^Xi, x, il.llb, ■ 

3 

All the maps are compatible with the respective symplectic, pre-sym- 
plectic and Poisson structures. Namely the map ^ is a composition of 
the quotient by the kernel of the pre-symplectic form and a symplectic 
map to a symplectic leaf. The map (|6| maps the symplectic form to the 
pre-symplectic one. The map ^ is Poisson, the map Q is anti-Poisson 
(Poisson with the opposite Poisson structure on the A'-torus). The maps 
([7| and ^ are dual to each other in the sense on Poisson pairs. 

Let I = (I, Io, e, d) and I' = (/', I' , e' , dl) be two seeds, and k 6 I — 1 . 
A mutation in the vertex k is an isomorphism fik '■ I — » I' satisfying the 
following conditions: 

• jti*(io) = ^0 = 

!— £ij if i = k or j = k otherwise 

£ij if e ik e kj < 

Eij + £ik\£kj\ if £ik£kj > 
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Two seeds related by a sequence of mutations are called equivalent. 
Mutations induce rational maps between the corresponding seed tori, 
which arc denoted by the same symbol /ik and are given by the formulae 



for the ^-torus, 



x k if i = k 

Xi(l + x k ) £ik if Sik > . 

Xi (l + {xk)- 1 )^ if e lk <Q 



for the ^4-torus and 



ll «T - n 

j\e jk >0 j|gjfc<0 

a-k 
a; 



if i = k 
iii^k 



{ (l+^r 1 n &r + (i+(x fe )- 1 )- 1 n &• 



i|ejfc>0 



j\ejk<0 



Vk 1 if i = k 

y l (l + y k ) £ ' k if e lk >0 

ViQ- + (Vk)- 1 )^ if e lk < 



if i = k 
iii^k 



for the 2?-torus. 

Since in the sequel we shall extensively use compositions of mutations 
called also cluster transformations we would like to introduce a shorthand 
notation for them. Namely, we denote an expression /U^.^)/^ by /Jj/x/s, 
M^. w)([ti (fe)M w (j)W by HkfijVi, and so on. 

Mutations have the following properties (valid for mutation of seeds 
as well as for mutations of respective tori) : 

• Every seed I = (I, I , e, d) seed is related to other seeds by exactly 
tt(7 — Jo) mutations. 
A\: Hi^i = id 
A\ x A\ If Eij = Sji = then fiiUj/ij/ii — id. 

A2: If Sij = —Sji = —1 then /i,/Xj/jj/ij/x, = id. (This is called the 
pentagon relation.) 



B2: If £%j 



G 2 : Ife 



*3 



-2eji = —2 then /Xi/ij/Xi/Xj/Xi/ij = id. 
-2>£ji = —3 then /jj/Xj/Xj/jj/ij/Xj/Ji = id. 
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By id we mean here an isomorphism of the seeds or tori. Conjecturally all 
relations between mutation follow from these ones. 

Given a seed one can produce a (t(7 — Iq) seeds by mutations. Contin- 
uing this procedure one obtains a §(I — Io)-valent tree whose vertices are 
seeds (or seed tori) and edges are pairs of mutually inverse mutations. Ob- 
viously if we start from any other seed from the tree we obtain the same 
tree. Every two tori of the tree are related by exactly one composition 
of mutations. Call two points of two different tori equivalent if they are 
related by the composition of mutations. The cluster manifold (denoted 
by A\i\ or T>\j\ depending on which kind of tori are used) is the affine 
closure of disjoint union of the tori quotiented by the equivalence relation. 

Each particular seed tori can be considered as a coordinate chart of 
the corresponding cluster manifolds and compositions of mutations can be 
considered as transition functions between the charts. 

Mutations respect the Poisson structure when acting on X tori, pre- 
symplectic structure when acting on ,4-tori and symplectic when acting on 
2?-tori. Thus the cluster manifolds X\i\, A\\\ and T>\\\ acquire the respec- 
tive structures. (In fact the formula for mutation of the matrix e can be 
considered as a corollary of this property and the mutation formulae for, 
say, X-tori). 

Mutations commute with the maps (|5| , ^ , ^ and ^ thus these maps 
are defined between the respective cluster varieties compatible with pre- 
symplectic, symplectic and Poisson structures thereof. 

Mutations are rational maps with positive integral coefficients and 
thus the cluster manifold can be defined not only over a field but over any 
scmificld as well. For scmificlds without -1 (like the semifields of positive 
real numbers or the tropical scmificlds) the mutations are isomorphisms 
and thus the whole manifold is isomorphic to every coordinate torus. 

The symmetry group D\i\ of a cluster manifold permuting the seed tori 
is called the (generalised) mapping class group of the cluster manifold. The 
name comes from the case of Teichmiiller space, when this group is the ac- 
tual mapping class group. The group depends on the equivalence class of a 
seed only and is common for cluster manifolds of types X , A and T>. Every 
sequence of mutations together with an isomorphism of the initial and the 
final seed gives an element of the mapping class group. Conversely, given 
a seed, every mapping class group element can be presented by a sequence 
of mutations starting from the given seed together with the isomorphism 
between the final seed and the initial one. Two sequences of mutations 
different by the relations A\-Gi correspond to the same mapping class 
group elements. 
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Consider the ring of algebraic functions on a cluster manifold in more 
details. The ring of algebraic functions 0(X{) (resp. 0(T>i), 0(Ai)) on 
every torus is the ring of Laurent polynomials of cluster variables. This 
ring contains a subring of Laurent polynomials with integral coefficients 
O z and a semiring of Laurent polynomials with positive integral coeffi- 
cients 0> o also depending of course of the seed and of the type of the 
torus. A cluster transformation in general does not presereve the ring O 
since it is birational. The ring of algebraic functions on the whole clus- 
ter manifold is the intersection of inverse images of the rings O under all 
possible cluster transformations of a seed tori. In other words the ring O 
consists of Laurent polynomials which stay Laurent under all possible clus- 
ter transformations. The celebrated result of Fomin and Zelevinsky called 
Laurent phenomenon claims that for any cluster variety A\\\ of type A the 
coordinate functions belong to the ring O 1 . The ring O contains a subring 
O z and a subsemiring 0> o - The latter is additively generated by Laurent 
polynomials from O with positive integral coefficients indecomposable into 
a sum of two such polynomials. Such elements of the semiring C>q £ir<3 
called irreducibles. The main conjecture, proven for a sufficiently wide 
class of cluster manifolds describes the structure of the set of irreducible 
Laurent polynomials: 

• The set of irreducible Laurent polynomials is a basis in the ring 
O. 

• The set of irreducible Laurent polynomials for the cluster variety 

(resp. ^4|i|, is canonically isomorphic to the set of points 
of the cluster variety A\i\{J}) (resp. ^(Z*), X>|i|(Z*)). 

One can consider this property as a duality between cluster varieties of 
type X (resp. A, V) and the tropical cluster varieties of type A, X and 
V, respectively. 

The correspondence between irreducible Laurent polynomials and 
points of the dual tropical variety is especially simple for the variety of 
type X. In this case the coordinates of the corresponding point of the 
tropical variety are given by multidegree of the highest term of the corre- 
sponding Laurent polynomial. 

Example 

Let us consider the simplest nontrivial example: the seed I — {I, e} 
with / = {1,2} and e 12 = 1. There are exactly 5 isomorphism classes of 
seed tori equivalent to a given one, however all the five seeds are isomor- 
phic, thus the mapping class group is Z/5Z. 
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The simplest geometric meaning has the space X. It is the space of 
5-tuples of points (p\, . . . ,p§) on the projective line P 1 such that pi ^ 
Pi+i (mod 5) an d modulo the automorphisms of P 1 . The 5-tuple of coordi- 
nate systems on this space is numerated by triangulations of the pentagon 
with vertices 1, . . . , 5. For every internal diagonal one associates the cross- 
ratio of the four points of the quadrilateral which this diagonal cuts into 
halves. Mutations correspond to removing a diagonal and replacing it by 
another one of the quadrilateral. The same variety over R>o is the config- 
uration space of 5-tuples of points on RP 1 with prescribed cyclic order. 

The „4-space is the space of collections of 10 nonvanishing vectors 
Vi,... ,vio in F 2 equipped with a nonzero bivector Vol. The collections 
are considered up to the action of the group SL(2,¥) of linear transfor- 
mations preserving Vol and subject to the relations = — v i+5 ( mod 10 ) 
and Vi A ( mo( j 10 ) = Vol. The map Xm — > Xm is given by the obvious 
projection of F 2 — {0} — > P 1 . For the internal diagonal of the pentagon 
with ends i and j one associates the coordinate v t Avj /Vol. The V variety 
is the space of flat SX(2,F) connections on a sphere with 5 different points 
on the equator removed with parabolic monodromy around these points. 
Consider the associated vector bundle and choose a monodromy invariant 
section about each singular points. Then trivialise the bundle over the 
northern hemisphere. The five chosen sections give five vectors vi, ... ,1)5 
in F 2 . The same procedure over the southern hemisphere gives five vectors 
W\ , . . . , w§ in another copy of F 2 . Given a triangulation of the pentagon 
we associate to every internal diagonal two coordinates x and b. The coor- 
dinate x is just the cross ratio of four points in P 1 defined by the vectors 
Vi standing at the corners of the quadrilateral cut by the diagonal (just 
like for the A'-space). The coordinate b is given by b — (vt A Vj)/(wi Awj), 
where i and j are the ends of our diagonal. The two projections to the 
X variety are obviously given by projectivising the collections of vectors 
{vi} and {wi}, respectively. The same manifold over IR>o can be identified 
with the space of complex structures on a sphere with five punctures on 
the equator. 

Given a triangulation of the pentagon one can describe the basis of 
the ring O z of the corresponding ^-variety explicitly as a set of Laurent 
polynomials P a ^{x,y) of two variables x, y parameterised by two integers 
a, b as follows: 
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x y" 

x 3 y b (l + x- 1 )- b 
x a y b (l + ar 1 )" b (l 



-l 



- x 1 y 1 ) a 

-l„,-lNa-b 



V 1 ) 



if 
if 
if 
if 
if 



a < 0, b < 
a < 0,b > 
a > 0,b < 
a > b > 
b > a > 



x a y b (l + y~ 1 )(l+ u 
x a y b (l + y~ 1 ) a 

One can easily check that this set of Laurent polynomials is invariant under 
simultaneous mutation of the variables x, y and of the variables a, b. 
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From max-plus algebra to non-linear Perron-Frobenius 
theory: an approach to zero-sum repeated gamea^ 

Stephane Gaubert 

This talk is based essentially on two joint works, with Akian and 
Nussbaum AGN07], on the one hand, and with Akian and Lem- 
mcns [AGLQ7], on the other hand. 

1. Introduction 

The analysis of zero-sum repeated games by the dynamic programming 
method classically leads to studying discrete time dynamical systems of 
the form 

(i.i) «(fc, .) = /(«(*- V)) 

where the map / is order preserving. Here, v(k, •) is the value function, 
which associates to any initial state the value of the corresponding game 
in horizon k. The map / is the "one day" dynamic programming operator. 
The case of a finite state space is already interesting. Then, denoting by 
n the number of states, we may identify the value function to a vector in 
K n , and the map / to a self-map of R". 



1 This work was partially supported by the joint RFBR/CNRS grant 05-01-02807 
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The explicit form of / depends on the details of the game. However, 
the map / may be written abstractly as: 

(1.2) f(v) = inf sup-/ 



where the infimum is taken over the strategics a of the first player and the 
supremum is taken over the strategies 7r of the second player, r aTC E 1" is 
a vector of payments, and P a7T is a n x n nonnegative matrix. In the case 
of games with undiscounted payoff, the matrices P a7T are stochastic. When 
there is a positive discount rate, or when the game may halt with a positive 
probability, the matrices P a ^ are substochastic. These (sub)-stochasticity 
properties imply that / is nonexpansive in the sup-norm, meaning that 

||/(«) - f(w)\\oo < \\v - w\\oo . 

The relevance of the order and nonexpansivencss properties to control 
and game problems has been brought to light by several authors, see in 
particular |CT80l IKol92L IRSOll |Ney03l 



The dynamic programming operators ( 1.2[ ) may be thought of as gen- 
eralizations of linear positive maps in several different ways. First, linear 
maps of the form x i— > Px, where P is a (sub)stochastic matrix, correspond 
to the zero-player case, in which every player has only one possible strat- 
egy, if we assume in addition that the payments are zero. Another special 
situation concerns the deterministic one player case, in which one of the 
two players has only one possible strategy, and the entries of the matrices 
P a7r are only or 1. Then, / becomes an affinc map over the min-plus 
or max-plus semiring. Further connections with Perron-Frobenius theory 
become apparent when using a familiar tropical instrument, the "loga- 
rithmic/exponential" glasses or "dequantization" , as in [LMSOll IVirOlj . 
This leads us to consider the conjugate map: 

g = expo/ o log- 
where log denotes the map from the interior of the standard positive cone 
R5 := {x € K™ | x > 0} to K n which does log entrywise, and exp := log -1 . 
Then, the map g is an order preserving self-map of the interior of K™ , and it 
is positively homogeneous or subhomogeneous of degree one, meaning that 
g(tx) = tg(x) or g(tx) < tg(x) for all scalars t > 1 and for all x e intM™. 
Such maps belong to non-linear Perron-Frobenius theory, which deals with 
the nonlinear extensions of the spectral theory of positive linear maps. We 
refer the reader to [Nus88 for a general account of this topic and for 
references. 

I will present some results concerning zero-sum games, which have 
been obtained by exploiting methods from non-linear Perron-Frobenius 
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theory with a max-plus or tropical point of view. The main results are 
taken from the two joint works |AGN07l [AGL07 . 

2. Nonlinear spectral radius of dynamic programming operators 

The classical notion of spectral radius has been extended to nonlinear maps 
in several ways |MPN02j . We assume here that g is a continuous posi- 
tively homogeneous of degree one map leaving invariant a (closed, convex, 
pointed) cone C in a Banach space X, and that g preserves the order in- 
duced by C, which is such that x < y if y — x € C. Bonsall's cone spectral 
radius of g is defined by: 

r c (g)=Mm\\g k \\y k 

where, for all continuous, positively homogeneous of degree one self-maps 
h of C, 

,,,,, I1M*)I1 

\\h\\c := sup . 

x£C\{0} \\ x \\ 

Another natural definition of the spectral radius arises when considering 
the nonlinear cigcnproblem: 

g{u) = \u 

where the nonlinear eigenvector u belongs to C \ {0}, and the nonlinear 
eigenvalue A is a nonnegative number. The cone eigenvalue spectral radius, 
fc(g), is by definition the maximal nonlinear eigenvalue A. Under some 
assumptions involving measures of non-compactness, it has been shown 
in [M PN02] that fc(g) = ?c(g)- Other useful notions of spectral radius, 
which coincide with the previous ones under reasonable assumptions, are 
studied in |MPN02j . 

We shall discuss here the related notion of Collatz- Wielandt number, 
which is obtained by considering super-eigenvectors in the interior of the 
cone instead of eigenvectors in the closed cone: 

?c{g) ■= inf{A > | 3u e int C, g(u) < \u} . 

The term "Collatz- Wielandt number" arises from Wielandt 's proof of the 
finite dimensional Perron- Frobenius theorem, in which the same formula is 
seen to characterize the Perron root of an irreducible nonnegative matrix. 

The main result of [AGN07 shows that rc(g) = r~c(g), when the 
cone C is normal, and when g satisfies some compactness assumptions. 



We apply these tools to dynamic programming maps of the form (1.2|, 
when the payments r a7r are 0, so that / is positively homogeneous of 
degree one. Under some standard assumptions (compactness of the action 
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spaces, continuous dependence of the reward and transition probabilities 
in the actions), which imply that the infimum and supremum are attained 



in (1.2| for all v, it is shown in |AGN07] that 



(2.1) f c (f) = r c (f) = r c (f) = inf sup r(P^) 

where C = M™ , and r(P cr7r ) denotes the Perron root of P CT7r . 

We derive from the previous result an explicit formula for the geomet- 
rical convergence rate of the iterates of the dynamic programming map /, 
this time with nonzero payments r a7T . To this end, we use the notion of 



subdifferential. Maps of the form ( 1.2 ) may not be differentiable, in partic- 



ular, if the action spaces are finite, / is piecewise affine. However, / may 
often be assumed to be scmidifferentiable, meaning that for all v and h, we 
can write f(v + h) = f(v) + f' v (h) + o(\\h\\), where f' v , the semidifferential 
of / at point v, is a continuous positively homogeneous of degree one map, 
which is defined uniquely by the latter property. 



When / is of the form ( 1.2 1, it can be shown that under fairly general 



assumptions, the semidifferential f' v at point v exists, and is given by: 
f v {h) = inf sup P a7T h , 



where denote the set of policies er which attain the infimum in ( 1.2 ) 



and for all a, H*(v, a) denotes the set of strategies tt which attain the 



supremum in the internal term in ( 1.2 1. (We need not assume that the inf 
and sup commute.) 

We show that if / has a fixed point v £ R™ , and if 

p := Tasx(rcU'v)i r -cU'v)) < 1 > 

then any orbit of / converges to v at a geometric rate which is bounded 
from above by p (this bound is tight). We eventually get the following 
explicit convergence rate: 

p = max( inf sup r(P CT7T ), sup inf r(P CT7r )) . 
3. Order preserving convex functions 

The techniques of the previous section are mostly appropriate when / has 
a unique fixed point, perhaps up to an additive or multiplicative constant. 
Therefore, a basic problem is to give a complete description of the 



fixed point set of the map ( 1.2 ). As a partial answer, a precise description 



of the set of stable fixed points is given in |AGL07] , when the map / is 
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convex (this corresponds to the one player case). This extends our ear- 
lier results [AG03J which concerned the undiscounted case. Here, we do 
not require any more the matrices P CT7r in (1.2 1 to be (sub) stochastic. In 
other words, we allow the possibility of a negative discount rate. Despite 
its apparently unphysical nature, negative discount is of practical interest: 
for instance, the study of static analysis problems by abstract interpreta- 
tion [GGTZ07 leads to fixed point problems involving maps which are 
always order preserving but not necessarily nonexpansive in some norm. 
Another motivation may come for fixed point problems for polynomials 
with positive coefficients, leading to maps like: 

fi(v) = log( ^2 a ij ex P0' ' 

where for all 1 < i < n, (aij)jeN" is an almost zero family of real nonneg- 
ative numbers. 

A convenient notion of stability, in the present setting, is the following 
one: we say that a fixed point v is *- stable if every orbit of the semidif- 
ferential f' v is bounded from above. It can be checked that a Lyapunov 
stable fixed point is ^-stable. 

Recall that a (communication) class of a nonnegative matrix P is by 
definition a strongly connected component of the digraph of P. We say 
that a class is critical if the corresponding principal submatrix of P has 
Perron root 1. The critical graph of P is the union of the subgraphs of the 
graph of P induced by the critical classes. If v is a ^-stable fixed point of 
/, we define the critical graph of /, G c (f) to be the union of the critical 
graphs of the matrices in the subdifferential 

df(v) := {P | f(w) - f(v) > P(w - v),Vw} . 

Of course, df(v) depends on v, but G c (f) is independent of the choice of 
the *-stable fixed point v. The critical nodes of / are defined to be the 
nodes of G c (/). 

We show in [AGL07] that a ^-stable fixed point is uniquely determined 
by its restriction to the set of critical nodes. Moreover, the restriction to 
the critical nodes allows us to identify the set of *-stable fixed points of / 
to a convex inf-subsemilattice of W, where p is bounded by the number of 
critical nodes. Some dynamical information, including a characterization 
of the possible lengths of "*-stable" periodic orbits of /, is also derived 
in |AGL07j . 

The results of |AG03l IAGL07] concern the one player case but have 
applications to the two player case. Indeed, the representation of the 
fixed point set has been used in [CTG06 to design a policy iteration 
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algorithm for zero-sum two player stochastic games. It allows one to 
handle "degenerate" iterations, in which the policies which arc selected 
yield dynamic programming maps with several fixed points. Some other 
applic ations o f these ideas, to static analysis of programs, are presented 
in |CGG+05LIGGTZ07| . 
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Cyclic projectors and separation theorems in 
idempotent semimodules^ 

Stephane Gaubert and Sergei Sergeev 
1. Introduction 

In an idempotent semiring, there is a canonical order relation, for which 
every element is "nonnegative" . Therefore, idempotent semimodulcs have 
much in common with the scmimodules over the semiring of nonnegative 
numbers, that is, with convex cones [6]. One of the first results based 
on this idea is the separation theorem for convex sets over "extremal alge- 
bras" proved by K. Zimmermann in [§]. Generalizations of this result were 
obtained in a work by S.N. Samborskh and G.B. Shpiz [7] and in works 
by G. Cohen, J. -P. Quadrat, I. Singer, and the first author pQ, [2]. 

The main result of this paper, Theorem |4.3[ shows that in the setting 
of finite-dimensional semimodules over max-plus semiring, several closed 
subsemimodules which do not have common nonzero points can be sepa- 
rated from each other. This means that for each of these subsemimodules, 
we can select an idempotent halfspace containing it, in such a way that 
these halfspaces also do not have common nonzero points. 

Even in the case of two semimodules, this statement has not been 
proved in the idempotent literature. Indeed, the earlier separation the- 
orems deal with the separation of a point from an (idempotent) convex 
set or semimodule, rather than with the separation of two convex sets or 
semimodules. 

In order to prove the main result, Theorem |4.3[ we investigate the 
spectral properties of idempotent cyclic projectors. By idempotent cyclic 
projectors we mean finite compositions of certain nonlinear projectors on 
idempotent semimodules. The continuity and homogeneity of these nonlin- 
ear projectors enables us to apply to their compositions, i.e. to the cyclic 
projectors, a result of R.D. Nussbaum[5] (non-linear Perron-Frobenius 

Supported by the RFBR grant 05-01-00824 and the joint RFBR/CNRS grant 
05-01-02807. 
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theory). We also show that the orbit of an eigenvector of a cyclic pro- 
jector maximizes a certain objective function. We call this maximum the 
Hilbert value of semimodules, as it is a natural generalization of Hilbert's 
projective metric, and characterize the spectrum of cyclic projectors in 
terms of these Hilbert values (Theorem 4.2 1. 

Our main results apply to the finite-dimensional semimodules over 
max-plus semiring. Some of our results still hold in a more general setting, 
see Sect. [3] However, the separation of several semimodules in such a 
generality remains an open question. 

The results of this paper are presented as follows. Sect. [2] describes 
the main assumptions, and some preliminary notions and facts that will 
be used in the paper. Sect. [3] is devoted to the results obtained in the 
most general setting, with respect to the assumptions of Sect. [2j The 
main results are obtained in Sect. [4] They include separation of several 
semimodules and characterization of the spectrum of cyclic projectors. 

The proofs of our results are contained in [3], which is an extended 
version of this text. 



2. Preliminaries 

We recall that a semiring (essentially, a ring without subtraction) is called 
idempotent, if its addition © is idempotent:a © a = a. The order relation 
mentioned above is given by a(Bb = b a < b. An example of idempotent 
semiring that will be important to us is K ma x.m, it is the set of nonnegative 
numbers M + equipped with operations a© b := max(o, b) and a© b = axb. 
It is isomorphic to the max-plus semiring (the set R U — oo equipped with 
a © b :— max(a, b) and a b = a + b). The "spaces" over semirings are 
called semimodules. 

An idempotent semiring or an idempotent semimodulc will be called 
b-complete, following [4], if it is closed under the sum (i.e. the suprcmum) 
of any subset bounded from above, and if the multiplication distributes 
over such sums. We shall consider semirings JC and semimodules V over JC 
that satisfy the following assumptions: 

(AO): the semiring JC is a 6-complete idempotent semifield, and the 
semimodule V is a 6-complete semimodule over JC; 

(Al): for all elements x and y ^ from V, the set {A S JC \ Xy < x} 
is bounded from above. 

Note that both assumptions are true for the semimodules JC 1 of re- 
valued functions on a set /, where JC is a fe-complctc semifield. 
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Assumptions (A0,A1) imply that the operation 

(2.1) x/y — max{A £ JC \ Xy < x}. 

is defined for all elements x and y ^ from V. 

Definition 2.1. A subsemimodulc V of V is a b-(sub)semimodule, if 
y is closed under the sum of any of its subsets bounded from above in V. 

Let V be a fo-subsemimodule of the semimodule V. Consider the op- 
erator Py defined by 

(2.2) P v (x) = max{ti e V \ u < x}, 

for every element x E V. Here we use "max" to indicate that the least 
upper bound belongs to the set. The operator Py is a projector onto the 
subsemimodule V, as Py(x) £ V for any x £ V and Pv(v) € V for any 

v e V. 

In idempotent geometry the role of halfspacc is played by the following 
object. 

Definition 2.2. A set H given by 

(2.3) H = {x \u/x >v/x}U {0} 

with u,v € Rmaxmi u — v i be called (idempotent) halfspace. 

Any halfspace is a semimodule. If V = JC n , an n-dimcnsional semi- 
module over /C, and all coordinates of u and v are nonzero, then we have 
that 

(2.4) H = {x\ x % u- x < xivr 1 }. 

{l,...,n} {1,...,"} 

The following theorem is a version of idempotent separation theorems 
[U[2], see also @]. 

Theorem 2.1. Le< V be a b-complete subsemimodule of V and let 
ueV be not in V. Then the set 

H = {x | P v (u)/x > u/x} U {0} 

contains V but not u. 

For any subsemimodule V and y G V, we denote 

(2.5) V v = {x & V \ y/x > 0}. 
It is a subsemimodule of V^. 
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Definition 2.3. A vector x is called archimedean, if x/y > for all 
y G V. A subsemimodulc of V is called archimedean, if it contains archi- 



medean vectors. A halfspace H defined by (2.3) will be called archimedean 
if both u and v are archimedean. 



Obviously Def. |2.3| makes sense only under 
(A2): The semimodulc V has an archimedean vector. 
This assumption is true in particular for the scmimodules of type /C™. 
In these semimodules we have that y/x > if and only if the support of 
x, i.e. the set supp(x) = {i | Xi ^ 0}, is a subset of supp(y) (the support 
of y). In this case V v has the form 

(2.6) V M = {xeV\ supp(x) C M}, 

for some index set M. A vector in JC n is archimedean if and only if it is 



positive. Regular halfspaces in this case are given by (2.4 1 



3. General results 

We shall study cyclic projectors, that is, compositions of projectors 

Pv k ---P Vl , 

where V\, . . . , V k are 6-subsemimodules of V. We assume (AO, Al), which 
means in particular that K, is an idempotent semifield. For the notational 
convenience, we will write Pt instead of Py t . We will also adopt a conven- 
tion of cyclic numbering of indices of projectors and semimodules, so that 
Pi +k = Pi and V i+k = Vi for all I. 

Definition 3.1. Let x 1 , . . . ,x k be nonzero elements of V. The value 

daix 1 , . . . , x k ) = (x^x 2 ) (x 2 /x 3 ) . . . ^/x 1 ). 

will be called the Hilbert value of x 1 , . . . , x k . 



The Hilbert value of two vectors x , x was studied in pQ. For two com- 
parable vectors in M^ axm , that is, for two vectors with common support 
M it is given by 

du{x\x 2 ) = mmixUx 2 )-^ 2 ^])- 1 ), 

so that — log(dYi(x 1 , x 2 )) coincides with Hilbert's projective metric 
5 n (x\x 2 )^\og(n^(x\(x 2 r 1 x 2 (x))- 1 )) = - \og{d u {x\ x 2 )). 
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Definition 3.2. The Hilbert value of k subsemimodules Vi, . . . , 14 of 
V is defined by 

d H (Vi, . . . , Vfe) = sup d B _(x 1 ,...,x k ) 

x 1 ev 1 ,...,x k ev k 

We establish two results on the spectrum of cyclic projectors and on 
their iterations. 

Theorem 3.1. Suppose that the operator P^o. . .oP\ has an eigenvector 
y with eigenvalue X, and define x l — Pi o . . . o P±y. Then 

\ = d H (V?,...,VV)=d H (x\...,x k ). 

Theorem 3.2. For any sequence of nonzero vectors {x l , i — 1, . . .} 
such that x 1 € V\ and x l = PiX 1 ^ 1 for i = 2, . . the Hilbert value 
dH(x l+1 , . . . , x l+k ) is nondecreasing with I. 



The following is an extension of Theorem 2.1 under assumptions (AO— 

A2). 

Theorem 3.3. Suppose that Vi, .. . , Vk are b-closed semimodules and 
that PfeO. ..oPj has an archimedean eigenvector y with nonzero eigenvalue 
X. The following are equivalent: 

(1) there exists an archimedean vector x and a scalar fi < 1 such 
that 

P fe o . . . o P x x < fix; 

(2) for all i = l,...,fc there exist regular half spaces Hi such that 
Vi C Hi and n.P, = {0}; 

(3) rw< = {0}; 

(4) A < 1. 



4. Projectors and separation in max algebra 

In M^ axm , it is natural to consider semimodules that are closed in the 
Euclidean topology. One can easily show that such semimodules are b- 
semimodules. Theorem 3.11 of [2] implies that projectors onto closed 
subsemimodules of M™ ax _ m are continuous. 

In order to relax the assumption concerning archimedean vectors in 
Theorem |3.3| we use some results from nonlinear spectral theory, that we 
next recall. By Brouwer's fixed point theorem, a continuous homogeneous 
operator x i— > Fx that maps R™ to itself has a nonzero eigenvector. This 
allows us to define the nonlinear spectral radius of F, 

(4.1) p(F) = max{A e R+ \ 3x € (R™) \ 0, Fx = Xx} . 
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Suppose in addition that F is isotone, then the maximum in (4.1l is at- 
tained and we can use the following nonlinear generalization of the Collatz- 
Wielandt formula for the spectral radius of a nonnegative matrix. 

Theorem 4.1. (R.D. Nussbaum, Theorem 3.1 of [5]) For any isotone, 
homogeneous, and continuous map F from Wt to itself, we have: 

p(F) = inf max \F{x)]iXj l . 

aG(R+\{0})" l<i<ri 

This result implies that the spectral radius of such operators is isotone: 
if F(x) < G{x) for any ielf, then p(F) < p{G). 

As the projectors on subsemimodules of K^ ax m are isotone, homogeneous 
and continuous, so are their compositions, i.e. cyclic projectors. Conse- 



quently, we can apply Theorem 4.1 to them. This allows us to refine 



the general results from the previous section. The following result refines 



Theorem 3.1 (the spectrum of cyclic projections). 



Theorem 4.2. Let Vi, . . . , Vk be closed semimodules in IJJj ax m . Then 
the Hilbert value of Vi, . . . , Vk is the spectral radius of Pfe o . . . o Pi . Ev- 
ery eigenvalue of Pk° ■ ■ ■ ° P\ is equal to d^V^ 1 , . . . , Vip ) for some M . 
Conversely, every such Hilbert value is an eigenvalue of Pk ° . . . ° Pi. 



The following result refines Theorem 3.3 (separation) 



Theorem 4.3. Suppose that Vi, i — 1, . . . ,k are closed semimodules 
o/M^ axm , and that fljVj = {0}. Then there exist archimedean half spaces 
Hi, i = 1, . . . , k such that Vi C Hi, for i = 1, . . . , k, and HiHi = {0}. 

A particular case of Theorem |4.3| is the following separation theorem 
for two semimodules. 

Theorem 4.4. Suppose that U and V are two closed max cones, and 
that U n V = 0. Then there exists an archimedean halfspace H\j , which 
contains U and does not intersect with V , and there exists an archimedean 
halfspace Hy , which contains V and does not intersect with U . 
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Pseudo-weak convergence of the random sets denned by a 
pseudo integral based on non-additive measure^ 

T. Grbic and E. Pap 



1. Introduction 

The weak convergence of sequence of probability measures is the main 
subject for a large class of limit theorems in the probability theory. In 
the classical probability theory it works with a-additive measures and 
the Lebesgue integral ([BJ. Several conditions equivalent to the weak 
convergence are provided by the theorem of Portmanteau (|B|). The main 
aim of this paper is to prove a Portmanteau-type theorem, with capacity 
functionals instead of probability measures, and with the general pseudo 
integral instead of the Lebesgue integral. 

Since the convergence in distribution of sequence of random closed sets 
on M can be tricky, it is often more appropriate to study the convergence 
of the corresponding sequence of capacity functionals. In this paper we 
study the convergence of sequences of random closed sets on R by looking 
at the convergence of the corresponding sequence of capacity functionals. 
Theoretical foundations of the theory of random sets, as generalization of 
random variables, were layed down by Kendall ([G]) and Matheron ([J]). 

Martially supported by the Project MNZZSS 144012, grant of MTA HTMT, 
French-Serbian project "Pavle Savic", and by the project "Mathematical Models for 
Decision Making under Uncertain Conditions and Their Applications" supported by 
Vojvodina Provincial Secretariat for Science and Technological Development. 
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Recall that random closed sets are random elements on the space of closed 
subsets of R. 

Our paper is organized as follows. Sect. 2 contains some preliminary 
notions, such as pseudo-operations and general pseudo integral [Al lOl 15]. 
In Sect. 3, we recall some basic notions and definitions from the theory 
of random sets ((Cj [Gj |Jl Kl |L]). The main results of this paper, also 
contained in Sect. 3, are concerned with the weak convergence of sequence 
of random closed sets, i.e., of the corresponding sequence of capacity func- 
tional with respect to the general pseudo integral. 



2. Preliminary notions 

Following [AjEllP]) we recall the notions of pseudo operations and general 
pseudo integral. Let < be the total order on [0,oo]. 

Definition 2.1. A binary operation © : [0,oo] 2 — > [0,oo] is called 
pseudo-addition if the following properties are satisfied: 
(Al) a © b = b © a ( commutativity) 
(A2) a<a'/\b<b'^a®b<a'e)b' (monotonicity) 
(A3) (a © b) © c = a © (6 © c) (associativity) 
(A4) a © = a (neutral element) 

(A5) a n — > a A b n — > b (a„ © 6„) — > a © 6 f continuity) 

Example 2.1. The following operations are pseudo-additions f [Al lOl 
IP]): (i) x © y = g^ 1 (g(x)) + g(y), where g : [0,oo] 2 — > [0, oo] is an 
increasing bijection; 

(ii) x © y — max(x, y) (note that this operation is idempotent). 

Definition 2.2. For a given pseudo-addition © pseudo-difference is 
the binary operation © : [0, oo] 2 — > [0, oo] given by 

a © b — ini{x <E [0, oo] : b © x > a}. 

Example 2.2. Obviously, a©6 = 0fora< & and a Q b > for 
a > b, see [Aj [I]. For pseudo-additions from Example 2.1 and a > b 
corresponding pseudo-differences are 

(i) aQb^g-^gia) - g(b)); (ii) aQb = a. 

Definition 2.3. For a given pseudo- addition © the pseudo-multipli- 
cation is a binary operation : [0, oo] 2 — > [0,oo] such that the following 
conditions are satisfied 

(Ml) a0O = O06=O (zero element) 
(M2) a < a' hb < b' ^ a Q b < a' Q b' (monotonicity) 
(M3) (a®b)Qc=(aQc)<3)(aQc) (right distributivity) 
(M4) o0l = l0d = a (unit element) 
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(M5) a Q (6 Q c) = (a &) c (associativity) 

(M6) a n — > a A b n — > (a„ 6 n ) — > a & f continuity) 



Example 2.3. (i) For the pseudo- addition from Example 2.1 (i), define 

pseudo-multiplication by a b = g~ 1 (g(a)g(b)). 

(ii) For the pseudo-addition from Example |2.1| (ii), one of the possible 

pseudo- multiplications is a 6 = a + b, see [E|, IF]. 

The algebraic structure ([0, oo],©,0) is a semiring. 

Let f2 be an abstract space, A a u-algebra of subsets of f2 and m : 
^t^Ra non-decreasing set function with m(0) = 0. We consider the space 
(tl,A,m) and a family of A-measurable functions / : f2 — > [0,oo], denoted 
by J 7 . A simple function is a measurable function s : — > [0, oo] whose 
range is finite. Let Rang(s) = {ai, ci2, . . . , a^} such that < a\ < a 2 < 
. . . < afe, and n Aj — for i 7^ j. The standard (B-step representation 

k 

of a simple function s is given by s — © b(c*,C*), where c* = a%, c* 2 — 

i=i 

m 

ffl2Q«i, . . . , c* m = a m a m _i, C* = ]J Aj and b : fl — > [0, 00] is a &asic 

i=« 



function of the form 6(c*, C*)(u>) 



0, g C? 



Definition 2.4. (i) The general pseudo integral of a simple function 
s with the standard ©-step representation is given by 

m m 

/ s0 dm= ©c* 0m(C*). 
i=i 

(ii) The general pseudo integral of a measurable function / G T is 
given by 

/* / © <^ m = SU P{/* 3 © ^ m : s G <5/}> 
where 5/ is the family of all simple function s such that s < f. 

The general pseudo integral has the following properties: 

(i) f® &(c, C) © dra=c0m(C). 

(ii) / <9=> /® / © rfm</%0 dm. 

(iii) For the pseudo characteristic function \A '■ fi — * [0,oo] of a set 

f 1 a; G A 

A C ft, defined by Xa{%) = \ n ' j 4 ' we have 



0, a; G 1 A 
J® dm = m(A). 
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3. Weak convergence of the sequence of capacity functionals 

3.1. Random closed sets and capacity functionals. We start with 
a short overview of the theory of random closed sets ([Cl |Gl |Jl |Kj iLl 
IMl IN]). Denote collections of closed, open and compact subsets of R by 
T ', O and /C, respectively. A very important role in the theory of random 
closed sets is played by collections of closed sets T, and its sub-collections 
T G = {F G T : F n G ^ 0}, G G O, and T K = {F G T : F n K = 
0}, K G K. Collections {T G : G G O} and {T K : K <E K,} generate a 
topology r(^ r ) on T . This topology is known as hit-or-miss-topology. The 
collection T endowed with the hit-or-miss topology is a compact, separable 
and Hausdorff space (jJ]). Taking countable unions and intersections of 
open sets of the topological space (T, t(J 7 )), we obtain a cr-field E(.F). 

Definition 3.1. A random closed set S is a measurable mapping from 
the probability space (SI, ^4, P) into the measurable space (J 7 , £(.F)). 

A random closed set S generates a probability distribution Ps in the 
following way 

P S (A) = P({lu € : S(w) € A}) = P S (S e A), for all A e E(^). 

Definition 3.2. For a random closed set S its capacity functional 
T s : K -> [0, 1] for X e /C is defined by 

T S (K) = P S (S G ^jf) = P S (S n K + 0). 

The capacity functional Tg is defined on /C, and it can be extended 
onto the family V of all subsets of K. A subset M C R is called capacitable 
if the following equality T S (M) = sup{T s (i ; s:) : K € K,, K C M} is true. 
All Borcl sets i? are capacitable ([L]). For a given random closed set S, and 
a sequence of random closed sets {S„} corresponding capacity functionals 
will be denoted by T and {T n }, respectively. 

3.2. (©, ©)-weak convergence. 

Definition 3.3. A sequence of capacity functionals {T„} (©, Q)-weak 
converges to a capacity functional T (shortly, pseudo-weak converges) if and 
only if for each continuous, bounded function / : R — > [0, oo] we have that 
lim f f dT n = f f dT. 

n — >oo 

We have proved in [D the following three theorems. 

Theorem 3.1. If a sequence of capacity functionals {T n } pseudo-weak 
converges to capacity functional T, then limsup T n (F) < T(F) for all 

n 

closed sets F CM.. 
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Theorem 3.2. If a sequence of capacity junctionals {T n } pseudo-weak 
converges to capacity functional T, then liminfT„(G) > T(G) for all 

n 

open sets Gel. 

Theorem 3.3. If for a sequence of capacity functionals {T„} and for 
all closed sets F holds (A) limsupT n (F) < T(F) and for all open sets 

n 

G holds (B) lim inf T n (G) > T(G), then {T n } pseudo-weak converges to 

n 

capacity functional T. 

Corollary 3.1. For a random closed set S and a sequence of random 
closed sets {S„}, which are defined in the following way: S = {X} and 
S n = {X„} ; where X is a random variable and {X„} is a sequence of 
random variables, the (®,Q)-weak convergence is equivalent to the weak 
convergence (with respect to continuous, bounded function f : R — > [0,00]/ 

Proof. For S = {X} and S„ = {X„}, we have that T(K) = P(X £ K) 
and T n (K) = P(X„ £ K) ([C]), where T and T„ are capacity functionals 
of random sets S and S„ , respectively. Since for each Borel set B we have 
that T(B) = sup{T(iT) : K £ /C, K C B}, it follows that T(B) = P(X € 
B). For all n e N we have that T n (B) = P(X„ e B). Suppose that {S„} 
(©,0)-weak converges to S. Then by Theorem 



3.1 



limsupP(X„ ef)< 



P(X € F) for all closed sets F. From the classical theorem of Portmanteau 
(|E]) we obtain that J fdP n — > J fdP, i.e., that the sequence of random 
closed sets {S„} weak converges to S. 

The weak convergence of the sequence of probability measures for any open 
set G implies that lim inf P(X„ e G) > P(X £ G), and for any closed set 

n 

F it implies that lim inf P(X„ £ F) < P(X £ F). Then, by Theorem 

n 

the sequence of random closed sets {S n } (®,©)-weak converges to S. 



3.3 



Remark 3.1. (i) For the special case described by corollary 3.1 the 
capacity functional reduces to the probability measure and then Theorem 
3.3| can be proved by taking into the account only one of the assumptions, 
(A) or (B). 

(ii) Weak convergence of the sequence of capacity functionals with 
respect to Choquet integral is investigated in [M_. Some equivalent con- 
ditions for the weak convergence of the sequence of probability measures, 
induced by sequence of random closed sets, are obtained in (Nl [Q] . 

(iii) The results obtained in this paper will serve for the investigation 
of further convergence properties of a sequence of capacity functionals of 
the sequence of random closed sets, based on the idempotent sup-measure 
and related integrals, see ([Dj |R] , 
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The stationary phase method and large deviations 

Oleg V. Gulinsky 

Let {Pa} be a family of probability measures on a measurable space 
(X, and let J be a nonnegative function on X with compact level sets. 
{P\} obeys the large deviation principle with a rate function / if and only 
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if 

Urn [f (g(x)) x P x (dx)] 1/X = S upg(x)e- 1 , 

A ^°° JX x£X 

for all bounded continuous nonnegative functions g on X pQ, [4] , 

We say that in this sense {Pa} converges to an idempotent measure 

exp{— /}. The r.h.s. of the last display is called a sup - integral or 

idempotent integral with respect to the idempotent measure and defines 

rough logarithmic asymptotics of the Laplace method. 

In this report we discuss logarithmic asymptotics of the integral 

J(A) = / exp{i\u(x)}(g(x)) x P x (dx), 
J x 

where X = R d (in what follows R 1 for the simplicity) and u,g (g > 0) are 
smooth enough real- valued functions. 

We consider this problem as a natural generalization of the station- 
ary phase method which imbeds the classical one in the context of large 
deviations. The interest in the problem is motivated by the slicing approx- 
imation approach to infinite dimension oscillatory integrals as well (see, for 
example [2]). 

Our approach is based on the technique of an almost analytic extension 
and follows the ideas of [3j where the classical method of stationary phase 
was extended to the case of complex-valued phase function. The new 
difficulty in our problem is the following. The function which plays the 
role corresponding to the imaginary part of the phase function in 3 , is just 
the rate function / defined asymptotically by the large deviation principle. 

Nevertheless, we consider f(x) — u(x) + il(x) as complex-valued 
"phase function" and assume that f(x) is C°° function in a neighbor- 
hood of the origin, which in turn is a non-degenerate stationary point of 
/ with 1(0) = 0. 

We introduce an almost analytic extension of / as follows: 

f(z) = [u(x)x(y) - l'(x)y X (t 1 y) - ^u" (x)y 2 x(t 2 y) + ■■■} 

+i[I{x) X {y) + u\x)yx{hy) - ^I"{x)y 2 x{t 2 y) + ...] = 

u(x,y) + iv(x,y), 

where x{y) G Co° IS ec l ua l to one in a neighborhood of the origin and 
vanishes for | y |> 1. The numbers tk > 1 are chosen sufficiently large so 
that the series converges. 
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To examine the asymptotic behavior of the integral J(A), we replace 
the integration along R by the integration along a suitable chain in the 
complex domain passing through the critical point of f(z). We show that 
on this chain the problem is reduced to the standard variational principle 
of large deviations. 

To fulfil the program following [3 , we first find new coordinates z in C 
for which f(z) — /(0) is a quadratic form in z . To this end, using Taylor's 
formula we write 

/(*)-/(0) =< z,R(z)z>/2, 

where by definition R(z) is an almost analytic function of z. Since all non- 
dcgcncratc quadratic: forms on C are equivalent, there is a linear transfor- 
mation A such that 

A T R(0)A = il. 

In turn, the equation 

iQ T (z)Q T (z) = R(z), 
Q(0) = A -1 , 

has a C°° solution Q(z) defined near the origin, since the map Q — > iQ T Q 
is analytic with surjective differential at Q — A -1 . Moreover, Q is an 
analytic function of R(z) and therefore an almost analytic function of z. 

The map z — > z = z(z) — Q(z)z defines new coordinates in the neigh- 
borhood of the origin and in this coordinates we have 

f(z) =f(0)+i<z,z>/2, 

where z — x + ly and < z, z >— x 2 — y 2 + 2i < x, y >. 

Since I(x) > with 1(0) = 0, it follows that x 2 ~ y 2 > on the 
tangent space at the origin and so there is a C°° function (p, defined in 
a neighborhood of 0, such that in the new coordinates R is given by the 
equation y = ip(x). 

We are now in a position to define a family of chains T s in C and 
examine the behavior of / on them. Let z — > z = z(z) be the inverse of 
the map z — > z(z). For < s < 1, putting 

T s = {z : z = z(z s ), z s — z s (x) — x + itp(x)s, x G R} 

one gets the estimate 



Imf(z(z s )) > C | Imz(z s ) \ 2 
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To replace the integration, we first consider the chain Ti and note that in 
a small enough neighborhood of z = the integrals J(A) and 

f exp{i\f(z)}(g(z)) x exp{XI(x) X (y)}I\(dx)dy, 

where g(z) is almost analytic extension of g, are equivalent( we may assume 
that the support of g w.r.t. z belongs to a small fixed neighborhood of the 
origin) . 

Finally we have to show that J Ti differs from J r with a very small 
error. We are able to do that with the help of Stokes's formula by the 
following arguments: (l)f(z) and g(z) are almost analytical functions, 
(2)Im/(*(*.)) > C | Imz(z s ) | 2 . 

Thus, it suffices compute the logarithmic asymptotics of the integral 

/ cxp{i\f(z)}(g(z)) x exp{XI(x) X (y)}Px(dx)dy = cxp{i\u(0)} x 

J exp{-A[/(0)+ | x | 2 /2}}(g(z(x)) x cxp{XI(x(x))}G(x)P x (dx), 

where G(x) = det(4i). One can easily recognize that the asymptotics of 
the last integral coincides with the asymptotics of 

exp{t\u(0)} x J (g(z(x)) x G(x)P x (dx). 

Thus we reduced the initial problem to the problem of large deviations. 
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Quantization with a deformed trace 

Dmitry Gurevich 

The standard quantization scheme of a Poisson structure on a variety 
M consists in the following. First, one looks for an associative ^-product 
satisfying the so-called correspondence principle. Existence of such a prod- 
uct is shown by Kontsevich. Second, one represents the constructed asso- 
ciative algebra An in a linear (hopefully, Hilbert) space. However, if the 
initial Poisson structure is not symplectic, such a representation is usually 
associated to each symplectic leaf of the bracket. 

In the 80 's the author considered some Poisson pencils whose quantiza- 
tion leads to "braided" algebras (cf. |Gll IG2] and the references therein). 
This means that in a sense they are related to a braiding, i.e. a solution 
to the Quantum Yang-Baxter equation (YBE) 

R12R23R12 = R23R12R23, where R 12 = R® I, R23 = I <S> R, 

V is a vector space over the ground field K (K or C), and R : V® 2 — > V® 2 
is a linear operator. Such a braiding plays the role of the usual flip in 
all related constructions and operations. In particular, generalized Lie 
algebras and their enveloping algebras were defined in this way. However, 
braidings entering their definitions were assumed to be involutive (R 2 = I) . 

Semiclassical counterpart of such braiding is a classical r-matrix. Given 
a Lie algebra g. By a classical r-matrix we mean an element r £ f\ 2 g sat- 
isfying the classical analog of the YBE 

[ru, ris] + i r i2, r 23 ] + [ri3, ^3} = 0, where r 12 =r®l,r 2 3 = l®r. 

Let p : g — * Vect (M) be a representation of the Lie algebra g into the 
vector fields space on a variety M. It is clear that the operator 

/®S6 K[M] -> {/, g} r = o p® 2 (r)(f <g> g) 

where o stands for the usual (commutative) product in the coordinate ring 
K[M] of the variety M defines a Poisson bracket on it. 

A typical example is M = g*, K[g*] = Sym(g). Given a classical 
r-matrix r € /\ (g) then the bracket { , } r is compatible with the linear 
Poisson-Lie bracket { , }pl, i-e. these brackets generate a Poisson pencil 

{, }a,6 = a{, }pL + b{, } r - 

Moreover, each of them (and consequently, the whole Poisson pencil) can 
be restricted to any G-orbit O C g*. (The restriction of the PL bracket 
to the orbit O is called Kirillov-Kostant-Souriau bracket.) 
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A quantization of the Poisson-Lie bracket can be realized in different 
ways. We consider the enveloping algebra U(gn) to be quantum coun- 
terpart of the bracket {, }pl- Hereafter by gn we mean the Lie algebra 
with the bracket K[ , ] where [ , ] is the bracket of the Lie algebra g and 
h is a deformation (quantization) parameter. As for the the KKS bracket 
its quantization can be realized as an appropriate quotient of the algebra 
U(gri) an d represented in a vector space in the spirit of the Kirillov orbit 
method. 

In order to quantize the whole pencil { , } a ,b or its restriction to an 
orbit O we apply the following result of Drinfcld. There exists an element 
F v € U{g)®U{g) such that 

F v = 1 + v r + . . . , 

F V {X + Y, Z) F V {X, Y) = F V (X, Y + Z) F U {Y, Z), 

and 

(e<8>l)F„ - (l®e)F„ = 1, 

where e is the counit in /7(g). 

By using this clement ("quantor" according to Lychagin's terminology) 
it is possible to quantize the above Poisson pencil and all other operators. 
Say, by equipping the algebra An = U (g/j) with a new product 

f*n,ug = *HP® 2 (F u )(f®g) 

where the representation p — ad is naturally extended to the algebra U{gf L ) 
and */j stands for the product in this algebra we get a new associative 
algebra denoted An,v 

The aforementioned braiding can be introduced via the clement F v . 
Namely, we put R = F^ 1 P ' F v where P is the usual flip. It is clear that 
R is involutive. Moreover, it is subject to the quantum YBE, i.e. it is a 
braiding. 

By means of the quantor F v the category of finite dimensional modules 
of the algebra An can be converted into that of ^4^-ones. This category 
is monoidal tensor rigid. Let us consider an object V of this category and 
the corresponding object End (V) = V <E) V* of internal endomorphisms. 
There exists a map 

Tr R : End (V) -» K 

which is a deformation of the usual trace and is morphism in this category. 
Moreover, it is R-symmetric, i.e. 

Tr R (X o Y) = Tv R o R{X <g> Y) 
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where o stands for the usual product in the algebra End (V). In a sense it 
looks like a super-trace for which the role of R is played by a super- flip. 

So, by quantizing the above Poisson pencil and by considering repre- 
sentations of the quantum algebra we are forced to replace the usual trace 
by its braided version. According to |G1] the linear term of the defor- 
mation of of the trace can be treated as a cocycle on the Lie algebra g. 
So, the deformation procedure itself can be regarded as a quantization of 
this cocycle. (Note that the involution operator A — > A* must be also 
deformed.) 

Recently it was understood what is an analog of the above algebra 
An,v corresponding to a non-involutive braiding (of Hecke type) and what 
is its semiclassical counterpart. 

Let R : V® 2 — > V® 2 be a Hecke symmetry, i.e. a braiding which meets 
the Hecke relation 

(ql - R)(q~ 1 I + R) = 0, q e K is generic. 

The algebra generated by the unit and elements lj, 1 < i,j < n = dim V 
subject to the equation 

RL\RL\ — L\RL\R = h^RLi — _Lii?), 

where L = (I?) is the matrix with entries l\ and L\ = L® 1 is call modified 
Reflection Equation Algebra (mREA). 

If the Hecke symmetry R comes from the quantum group U q {sl{n)) it 
is a one parameter deformation of the usual flip. In this case the mREA is 
two parameter deformation of the commutative algebra Sym (gl{n)) (which 
is a specialization of the An, q at h = 0, q — 1) . Its semiclassical counterpart 
is a Poisson pencil similar to that above but with the bracket { , } r defined 
in another way. Namely, it is an extension to the ambient vector space of 
the so-called Semenov-Tian-Shansky bracket defined on the group SL(n). 
Similarly to the pencil above the latter one can be also restricted to any 
GL(n)-orbit in gl(n)* . 

The algebra An, q possesses a braided bi-algebra structure and has 
the same category of finite dimensional representations as the quantum 
group Uq(sl(n)) has. However, in contrast with the above monoidal tensor 
category this one is quasitensor one. Nevertheless, an intrinsic trace Tr^ 
which is a categorical morphism and a deformation of the usual trace is 
well defined on any object End ( V) of internal cndomorphisms. (For simple 
objects V it is unique up to a factor.) 
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For instance, if V is the basic space then Tr# l\ = Sf . The defining 
relations of the algebra An, q can be rewritten as follows 

X ® Y - Q(X ®Y) = [X,Y], I,FeL = span 

where Q : L® 2 — > L® 2 is a braiding and [ , ] : L® 2 — > L is a " braided Lie 
bracket" . We would like to emphasize that the latter form of the mREA 
makes it more similar to an enveloping algebra. 
It is easy to check that 

Tr R {X oY) = Tr R o Q(X ® Y) 

where o stands for the usual product in the algebra End (V) (note that 
the space L can be naturally identified with End (V)). So, we can see that 
such a trace is Q-symmetric. 

A more detailed presentation of the topic can be found in the paper 
[GPS] . 

In my talk I shall exhibit the role of a deformed (quantum) trace in 
" braided geometry" . 
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Transformations preserving matrix invariants over 

semirings^ 

Alexander E. Guterman 

The investigations of matrix transformations which leave fixed differ- 
ent matrix properties and invariants is an actively developing part of ma- 
trix theory. This research was started in the works by Frobenius, see [gl 
Theorem 1.1], and Dieudonne, see [5j[9l Theorem 1.2], where bijective lin- 
ear transformations on matrices over fields which preserve the determinant 
and the set of singular matrices, correspondingly, were characterized. 

1 Partially supported by the RFBR grant 05-01-01048 and the grant MK- 
2718.2007.1. 



Transformations preserving matrix invariants over semirings 



85 



During the last three decades many authors investigated linear trans- 
formations on more general algebraic structures, such as matrices over 
rings and semirings. In this talk we are going to discuss the corresponding 
problems on max-algebras and related classes of semirings. 

Definition 4. A semiring is a set S with two binary operations, 
addition and multiplication, such that: 

• S is an abelian monoid under addition (identity denoted by 0); 

• S is a semigroup under multiplication (identity, if any, denoted 
by 1); 

• multiplication is distributive over addition on both sides; 

• sO = 0s = for all s € S. 

In this paper we will always assume that there is a multiplicative identity 
1 in S which is different from 0. 

Definition 5. A semiring S is called commutative if the multiplica- 
tion in S is commutative. 

Definition 6. A semiring S is called antinegative (or zero -sum- free) 
if a + b = implies that a = b = 0. 

This means that the zero element is the only element with an additive 
inverse 

Definition 7. We say that a semiring S has no zero divisors if from 
ab = in S it follows that either a = or b = 0. 

Definition 8. A semiring is called a max-algebra if the set S is an 
ordered group with the multiplication * and the order relation <, and 
operations in S are defined as follows: a + b = max{a, b}, a ■ b = a * b for 
any a, b G S. 

It is straightforward to see that max-algebra is antinegative. Also it 
does not contain zero divisors, moreover any non-zero element of a max- 
algebra has a multiplicative inverse. 

Let M. m>n {S) denote the set of to x n matrices with entries from the 
semiring S, A4 n (S) = Ad n>n (S). Under natural definitions of matrix addi- 
tion and multiplication A4 n (S) is obviously a semiring. Matrix theory over 
semirings has been an object of intensive study during the last decades, see 
for example the monograph [7 and references therein. The development of 
linear algebra over semirings certainly requires such an important matrix 
invariant as the determinant function. However it turns out that even over 
commutative semirings without zero divisors the classical determinant can 
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not be defined as over fields and commutative rings. The main problem lies 
in the fact that in semirings which are not rings not all elements possess 
an additive inverse. A natural replacement of the determinant function 
for matrices over commutative semirings is the bideterminant known for 
many years, see [7]- 

Definition 9. A bideterminant of a matrix A = [cn^] & M n (S) is 
the pair (||A|| + , ||A||~), where 

H-4|| + = 51 a l,cr(l) ' ' ' a n,a(n), \\A\\~ = X] a l,<x(l) " ' a n,cr(n) : 




here S n denotes the symmetric group on the set {1, . . . , n}, A n denotes its 
subgroup of even permutations. 

It is known that the bideterminant function possesses some natural 
properties. Namely it is invariant under transposition, and for any scalar 
a G S, (||ayl|| + ,||aA||-) = (a n ||A|| + , a n ||A||-). However, some basic 
properties of the determinant are no longer true for the bideterminant. 
For example, if A is invertible then \\A\\ + ^ \A\~ but the converse is not 
always true. 

Example 1. Let us consider A = E 1A +2Ei_ 2 +3E 2 ^+4:E2^ G M 2 (<S), 
where S — (Q+ , max, • ) , namely the set of non-negative rationals with the 
standard multiplication and the addition defined by a + b = max{a,6}. 
Then (||A|| + , P||") = (4,6) but A is not invertible. 

Note that the bideterminant is not multiplicative in general. However, 
some weaker versions of this property are true, in particular, 



\\AB\\+ + \\A\\ + \\B\\- + \\A\\~ \\B\\+ = \\AB\\- + \\A\\ + \\B\\+ + \\A\\- \\B\ 



We prove the following theorem which is a semiring analog of famous 
Frobenius theorem, see [6] , on linear transformations preserving the deter- 
minant of complex matrices. 

Theorem 1. [2] LetS be a commutative antinegative semiring without 
zero divisors and T : M. n (S) — > M. n (S) be a surjective linear transforma- 
tion. Then (||T(A)|| + , ||T(X)||-) = (||A|| + , p"|| - ) for all X G M„(S) if 
and only if there exists permutation matrices P, Q and invertible diagonal 
matrices D,E, satisfying (||PQ||+, ||PQ|| _ ) = (\\DE\\ + ,\\DE\\-) = (1,0), 
such that either T(X) = PDXEQ for all X G M n {S) or T(X) 
VI) X l\( ) for all X G M n {S). Here the matrices P, Q are defined uniquely 
and the matrices D, E are defined uniquely up to an invertible scalar factor. 
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Definition 10. We say that a transformation T : M n {S) — > M n {S) 
is standard if it is denned by T(X) = PDXEQ for all X G M n (S) or 
T(X) — PDX l EQ for all X G M n (S) for certain permutational matrices 
P, Q and diagonal matrices D, E. 

The following similar function is widely considered in combinatorial 
matrix theory: 

Definition 11. A permanent of a matrix A = [ctij] G A4 n (S) is 



Also the following polynomial is related to this function: 

Definition 12. The rook polynomial of a matrix A 6 M m ^ n (S) is 
-R^cc) = ^2 Pj x3 ' i where Pq = 1, p 3 is the sum of the permanents of all 

j x j submatrices of A. 

Linear transformations preserving the rook polynomial and permanent 
itself were characterized by Beasley and Pullman. 

Here we provide more general result, namely, we prove that in order to 
characterize a transformation, it is enough to know that it preserves any 
single coefficient of a rook polynomial, namely we prove the following: 



Theorem 2. [8] Let T : M n (S) — > M n (S) be a surjective linear trans- 
formation and j , 2 < j < n, be fixed. Then T preserves the j-th coefficient 
of the rook polynomial iff T is standard with pj(DE) = 1. 

The following notions of singularity arc in use while dealing with ma- 
trices over semirings, as usual, we separate left and right singularity. 

Definition 13. A matrix A e M m ,n{S) is said to be S-right singular 
if Ax = for some nonzero x G S n . A G Ai m ,n{S) is S-left singular if 
x*A = 0* for some nonzero x G <S TO . A matrix A G M. m ,n(S) is S-singular 
if A is either 5-left singular or 5-right singular. 

The next example shows that even over antinegative commutative 
semirings without zero divisors there exist matrices that are 5-left sin- 
gular and are not 5-right singular or vice versa. 

Example 2. Let (M, +,max) be a max-algebra, 




A = 





1 1 



,B = 



1 
1 



G A4 2 (K,+,max). 



We have that Ax. = forces x = since (M,+,max) is antinegative, but 
[1,0] A = [0, 0]. Similar x l B = forces x = while B[0, 1]* = [0, 0]*. 
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Definition 14. A matrix A e M. m , n (S) is S-nonsingular if A is not 
5-singular. 

Note that if S is commutative and A is an 5-singular square matrix 
then (||j4|| + , ||j4||~) = (0,0). However the following example shows that 
there are 5-nonsingular matrices with the bideterminant equal to (0,0). 

Example 3. Over any commutative antinegative semiring, 









1 


+ 








1 


1 


1 





= = 


1 


1 











1 










1 



We obtain the following analog of Dicudonnc theorem on singularity 
preservers, see [5], for matrices over semirings. 

Theorem 3. [2] Let S be an antinegative semiring without zero divi- 
sors and T : M m ,n(S) — > M m ,n(S) be a surjective linear operator. Then 
the following statements are equivalent 

(1) T preserves the set of S -singular matrices; 

(2) T preserves the set of S-nonsingular matrices; 

(3) There are permutational matrices P,Q £ A4 n (S) and a matrix 
B with all invertible entries such that either T{X) = P(X o B)Q 
for all X e M n (S) or T(X) = P(X o BfQ for all X S M n {S), 
here X o B is an Hadamard product, i.e., (X o B)ij = Xijbij. 

If the semiring S is also a subsemiring of an associative ring 1Z without 
zero divisors we can consider the following notion of singularity as well. 

Definition 15. We say that a matrix A e M m ,n{S) is IZ-right singu- 
lar if Ak = for some nonzero x € lZ n . A € M. m ,n(S) is IZ-left singular 
if x'A = for some nonzero x 6 lZ m . A is IZ-singular if A is either 7\L-left 
singular or 7?.-right singular, and IZ-nonsingular if it is not 7?.-singular. 

It is straightforward to see that if a semiring S is a subsemiring of a cer- 
tain ring 7Z then 7?.-right (left) singularity follows from 5-right (left) singu- 
larity. However the following example shows that there are iS-nonsingular 
matrices which are 7?.-singular. 

n 

Example 4. For any n the matrix J n — G -Mn(Z+) is Z-left 

and Z-right singular but Z + -nonsingular. 

Note that similarly to the situation over fields all non-square matrices 
are 7?.-singular, however as the above example shows they may not be 
5-singular. 
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An analog of Theorem[3]holds for transformations preserving ^-singularity. 
Corresponding transformations appear to be standard. 

Definition 16. Let S be a max-algebra (operations are denoted by 
max and +). A matrix A — [ay] G M n (S) is said to be tropically singular 
if the maximum in the expression for the permanent 

per (A) = max{a lff m + • . . + a nCT(n) } 
is achieved at least twice. 

It can be generalized to the case of an arbitrary antinegative semiring 
S in the following way: 

Definition 17. A matrix A = [ay] e M n (S) is said to be tropically 
singular if there exists a subset T € S n such that 

22 a l<r(l) ' ' ' a na(n) = ^ a l<j(l) ' ' ' a n(T („) . 
ff£T aeS n \T 

Our further results include the characterization of linear transforma- 
tions preserving these and related notions of singularity. Also we obtain 
several analogs of Markus and Moyls result on linear transformations pre- 
serving rank, see [9j Theorems 3.1, 3.2], for several well-known semiring 
rank functions. 

Bibliography 

[1] L. B. Beasley, N. J. Pullman, Term rank, permanent and rook polynomial preservers, 

Linear Algebra Appl. 90 (1987) 33-46. 
[2] L. B. Beasley, A. E. Guterman, S.-G. Lee, S.-Z. Song, Frobenius and Dieudonne 

theorems over semirings, Linear and Multilinear Algebra, 55, no. 1, (2007) 19-34. 
[3] P. Butkovic, Max-algebra: the linear algebra of combinatorics? Linear Algebra 

Appl. 367 (2003) 315-335. 
[4] M. Develin, F. Santos, B. Sturmfels, On the rank of a tropical matrix. In Discrete 

and Computational Geometry (E. Goodman, J. Pach and E. Welzl, eds.), MSRI 

Publications, Cambridge Univ. Press, 2005. 
[5] J. Dieudonne, Sur une generalisation du groupc orthogonal a quatre variables, Arch. 

Math. 1 (1949) 282-287. 
[6] G. Frobenius, Vber die Darstellung der endlichen Gruppen durch lineare Substitu- 

tionen, Sitzungsber., Preuss. Akad. Wiss (Berlin), Berlin (1897) 994-1015. 
[7] K. Glazek, A Guide to the Literature on Semirings and their Applications in Math- 
ematics and Information Sciences, Kluwer Academic Publishers, 2002. 
[8] A. E. Guterman, Transformations of non-negative integer matrices preserving the 

determinant, Uspehi Mat. Nauk 58, no. 6, (2003) 147-148. 
[9] C.-K. Li, S. Pierce, Linear preserver problems, Amer. Math. Monthly 108, no. 7, 

(2001) 591-605. 



90 



I. Itenberg, V. Kharlamov, E. Shustin 



Tropical geometry and enumeration of real rational 

curved 

/. Itenberg, V. Kharlamov, and E. Shustin 

The talk is devoted to applications of tropical geometry in enumera- 
tive (complex and real) algebraic geometry. We concentrate ourselves at 
enumeration of real rational curves interpolating fixed collections of real 
points in a real algebraic surface E, and more precisely, at the following 
question: given a real divisor D and a generic collection w o/ci(£) • D — 1 
real points in E, how many of the complex rational curves belonging to the 
linear system \D\ and passing through the points of w are real ? By ratio- 
nal curves we mean irreducible genus zero curves and their degenerations, 
so that they form in \D\ a projective subvariety S(T,,D); this subvariety 
is called the Severi variety. A curve on a real surface E is called real, if 
the curve is invariant under the involution c : E — > E defining the real 
structure of E. 

While, under mild conditions on E and D, the number of complex 
curves in question is the same for all generic collections w (it equals to the 
degree of S(T,,D)), it is no more the case for real curves (except few very 
particular situations) . 

J.-Y. Welschinger [5j [6] discovered a way to attribute weights ±1 
to the real solutions in question so that the number of real solutions 
counted with weights becomes independent of the choice of a generic col- 
lection of real points. As an immediate consequence, the absolute value 
of the Welschinger invariant Ws t o provides a lower bound on the number 
Rs,d(w) of real solutions: Rz y rj(w) > \W^,d\- 

In some cases (for example, in the case of toric Del Pezzo surfaces; 
recall that there are five toric Del Pezzo surfaces: the projective plane P 2 , 
the product P 1 x P 1 of projective lines, and P 2 blown up at k points in gen- 
eral position, where k = 1, 2 or 3) Welschinger invariants can be calculated 
using Mikhalkin's approach [3J |4] which deals with a corresponding count 
of tropical curves. In tropical geometry, complicated non- linear algebro- 
geometric objects are replaced by simpler piecewise-linear ones. For ex- 
ample, tropical plane curves are piecewise-linear graphs whose edges have 
rational slopes. Tropical curves can be seen as algebraic curves over the 
tropical semi-field (max, +). 

Using the tropical approach, we proved (see pQ) the logarithmic equiv- 
alence for the Welschinger and Gromov-Witten invariants of any toric Del 
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Pezzo surface equipped with its tautological real structure, i.e., the real 
structure which is provided by the toric structure. 

Theorem 1. (see [1]) Let £ be a toric Del Pezzo surface equipped 
with its tautological real structure, and D an ample divisor on S. The 
sequences logWs, n D an d logGWs.nD, n £ N, of the Welschinger in- 
variants and the corresponding Gromov- Witten invariants are asymptot- 
ically equivalent. More precisely, log Ws,nD = log GWs,nD + 0{n) and 
logGW^ nD = (ci (£) • D)-n log n + 0{n). 

We also defined (see [2]) a series of relative tropical Welschinger- type 
invariants of real toric surfaces. In the Del Pezzo case, these invariants 
can be seen as real tropical analogs of relative Gromov- Witten invariants, 
and are subject to recursive formulas of Caporaso-Harris type. 

In the present talk, we consider generic collections of real points on the 
projective plane blown up at 4 real points in general position and prove 
that the logarithmic equivalence of the Welschinger and Gromov- Witten 
invariants holds in this situation as well. 

Theorem 2. Let £ be the projective plane P 2 blown up at 4 real 
points in general position, and D an ample divisor on £. The sequences 
logV^s.nD and logGW^nD; n G N, of the Welschinger invariants and the 
corresponding Gromov- Witten invariants are asymptotically equivalent. 

The proof is based on a new version of the correspondence theorem, 
whose proof in turn uses an appropriate tropical Caporaso-Harris type for- 
mulas. In particular, we get recursive formulas that allow one to calculate 
Welschinger invariants of P 2 blown up at 4 real points in general position. 
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Abstract convexity and cone-vexing abstractions 

Semen S. Kutateladze 

This talk is devoted to some origins of abstract convexity and a few vexing 
limitations on the range of abstraction in convexity. Convexity is a rela- 
tively recent subject. Although the noble objects of Euclidean geometry 
are mostly convex, the abstract notion of a convex set appears only after 
the Cantor paradise was founded. The idea of convexity feeds generation, 
separation, calculus, and approximation. Generation appears as duality; 
separation, as optimality; calculus, as representation; and approximation, 
as stability. 

1. Generation 

Let E be a complete lattice E with the adjoint top T :— +00 and bottom 
_L := —00. Unless otherwise stated, Y is usually a Kantorovich space which 
is a Dedekind complete vector lattice in another terminology. Assume 
further that H is some subset of E which is by implication a (convex) 
cone in E, and so the bottom of E lies beyond H . A subset [/of H is 
convex relative to H or H- convex, in symbols U £ 23(ii, E), provided that 
U is the H-support set Up :— {h g H : h < p} of some element p of E. 

Alongside the if -convex sets we consider the so-called if-convcx ele- 
ments. An element p G E is H-convex provided that p = supf/^; i.e., p 
represents the supremum of the iJ-support set of p. The iJ-convex ele- 
ments comprise the cone which is denoted by ^(H, E). We may omit the 
references to H when H is clear from the context. It is worth noting that 
convex elements and sets are "glued together" by the Minkowski diality 
ip : p 1— » Up . This duality enables us to study convex elements and sets 
simultaneously. 

Since the classical results by Fcnchel pQ and Hormander [2J |3] it has 
been well known that the most convenient and conventional classes of 
convex functions and sets are ^{A{X), R x ) and %}(X' ,R X ). Here A is a 
locally convex space, X' is the dual of X, and A(X) is the space of afhnc 
functions on X (isomorphic with X'xl). 

In the first case the Minkowski duality is the mapping / epi(/*) 
where 

f*(y) := sap({y,x) - f(x)) 
xex 

is the Young-Fenchel transform of / or the conjugate function of / . In the 
second case we prefer to write down the inverse of the Minkowski duality 
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which sends U in 9J(X',K ) to the standard support function 

<y5 _1 (J7) : x i * sup(y,x). 

yeu 

As usual, (•, •) stands for the canonical pairing of X' and X. 

This idea of abstract convexity lies behind many current objects of 
analysis and geometry. Among them we list the "economical" sets with 
boundary points meeting the Pareto criterion: capacities, monotone semi- 
norms, various classes of functions convex in some generalized sense, for 
instance, the Bauer convexity in Choquet theory, etc. It is curious that 
there are ordered vector spaces consisting of the convex elements with re- 
spect to narrow cones with finite generators. Abstract convexity is traced 
and reflected, for instance, in [4]-[9]- 

2. Separation 

Consider cones K\ and K 2 in a topological vector space X and put x := 
{K\,K2). Given a pair x define the correspondence Q x from X 2 into X 
by the formula 

$x := {(k u k 2 ,x) £ X 3 : x = k t - k 2 e K % (t := 1,2)}. 

Clearly, <J>^ is a cone or, in other words, a conic correspondence. 

The pair x is nonoblate whenever <f> x is open at the zero. Since 
&x(V) — V H Ki — V n K 2 for every V C X, the nonoblateness of x 
means that 

xv := (V n k x - v n k 2 ) n (V n k 2 - v n K±) 

is a zero neighborhood for every zero neighborhood V d X. Since xV C 
V — V, the nonoblateness of x is equivalent to the fact that the system of 
sets {xV} serves as a filterbase of zero neighborhoods while V ranges over 
some base of the same filter. 

Let A„ : x 1— * (x, . . . ,x) be the embedding of X into the diagonal 
A n (X) of X n . A pair of cones x := (K\, K 2 ) is nonoblate if and only if 
A := (Ki x K 2 ,A 2 (X)) is nonoblate in X 2 . 

Cones K\ and K 2 constitute a nonoblate pair if and only if the conic 
correspondence <i> C X x X 2 defined as 

$ := {(h, xi, x 2 ) E X x X 2 : x, + h £ K % (t := 1, 2)} 

is open at the zero. Recall that a convex correspondence <i> from X into Y 
is open at the zero if and only if the Hormander transform ofXx$ and 
the cone A 2 (X) x {0} x K + constitute a nonoblate pair in X 2 x Y x R. 
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Cones K\ and K 2 in a topological vector space X are in general 
position provided that 

(1) the algebraic span of K\ and K 2 is some subspace Xq C X; i.e., 
X = K 1 -K 2 = K 2 -K 1 ; 

(2) the subspace Xq is complemented; i.e., there exists a continuous 
projection P : X — > X such that P(X) — Xq; 

(3) JCi and _K" 2 constitute a nonoblate pair in X . 
Let <7„ stand for the rearrangement of coordinates 

o- n ■ ((xi,yi), ■ ■ ■ ),(yi,---,yn)) 

which establishes an isomorphism between (X x Y) n and X™ x F™. 

Sublinear operators Pi, . . . , P n : X — > _EU{+c»} are in general position 
if so are the cones A n (X) x B n and cr„(epi(Pi) x • ■ ■ x epi(P„)). A similar 
terminology applies to convex operators. 

Given a cone K C X, put 

n E (K) := {T G _£f(X, E):Tk<0(ke K)}. 

We readily see that -^(-ftr) is a cone in J£(X, E). 

Theorem. Let K\, . . . , -ftT„ be cones in a topological vector space X 
and let E be a topological Kantorovich space. IfK%,..., K n are in general 
position then 

tte(Ki H • •• n K n ) = -KsiK-i) + ■■■ + ir E (K n ). 

This formula opens a way to various separation results. 

Sandwich Theorem. Let P, Q : X — > E 1 U {+00} be sublinear op- 
erators in general position. If P(x) + Q(x) > for all x £ X then there 
exists a continuous linear operator T : X — > E such that 

-Q(x) <Tx< P(x) (x G X). 

Many efforts were made to abstract these results to a more general 
algebraic setting and, primarily, to semigroups. The relevant separation 
results are collected in |10j . 

3. Calculus 

Consider a Kantorovich space E and an arbitrary nonempty set 21. Denote 
by loo (21, E) the set of all order bounded mappings from 21 into E; i.e., 
/ G ZooCSl, E) if and only if / : 21 — ► E and the set {/(a) : a G 21} is order 
bounded in E. It is easy to verify that Zoo (21, E) becomes a Kantorovich 
space if endowed with the coordinatewise algebraic operations and order. 
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The operator £21, e acting from ?oo(2l, E) into E by the rule 

e%, E ■ f -> sup{/(a) : a e 21} (/ e loo(%E)) 

is called the canonical sublinear operator given 21 and E. We often write 
£a instead of £21, e when it is clear from the context what Kantorovich 
space is meant. The notation e n is used when the cardinality of 21 equals 
n and we call the operator e n finitely- generated. 

Let X and E be ordered vector spaces. An operator p : X — > E is 
called increasing or isotonic if for all x±,X2 £ -X" from xi < X2 it follows 
that p(xi) < p(x2). An increasing linear operator is also called positive. 
As usual, the collection of all positive linear operators in the space L(X, E) 
of all linear operators is denoted by L + (X,E). Obviously, the positivity 
of a linear operator T amounts to the inclusion T(X + ) C E + , where 
X + := {x e X : x > 0} and E + := {e e E : e > 0} are the positive 
cones in X and E respectively. Observe that every canonical operator is 
increasing and sublinear, while every finitely-generated canonical operator 
is order continuous. 

Recall that dp := dp(0) = {Te L(X, E) : (Vx e X) Tx < p{x)} is the 
subdifferential at the zero or support set of a sublinear operator p. 

Consider a set 21 of linear operators acting from a vector space X 
into a Kantorovich space E. The set 21 is weakly order bounded if the set 
{ax : a e 21} is order bounded for every x e X. We denote by (21) x the 
mapping that assigns the element ax € E to each a S 21, i.e. (21) x : a 1— » 
ax. If 21 is weakly order bounded then (21) x e Zoo (21, E) for every fixed 
x e X. Consequently, we obtain the linear operator (21) : X — > ?oo(2l, -E 1 ) 
that acts as (21) : x (21) x. Associate with 21 one more operator 

P21 : x 1— » supjax : a G 21} (x G A). 

The operator P21 is sublinear. The support set dp% is denoted by cop(2l) 
and referred to as the support hull of 21. These definitions entail the 
following 

Theorem. If p is a sublinear operator with dp — cop(2l) then P = 
£21 o (21). Assume further that p\ : X — > E is a sublinear operator and 
P2 : E — > -F is an increasing sublinear operator. Then 

d(p 2 op 1 ) = {To(dp 1 ):TeL+{l 00 {dp 1 ,E),F) A T o A dpi e dp 2 } . 

Furthermore, if dp\ = cop(2li) and dpi = cop(2l2) then 

0(P2 o pi) = {T o (210 = T e i+(/oo(ai, £), F) 
A (3a e 9£2i 2 ) To A ai = ao (2t 2 )} • 
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More details on subdifferential calculus and applications to optimality 
are collected in 

4. Approximation 

Study of stability in abstract convexity is accomplished sometimes by intro- 
ducing various epsilons in appropriate places. One of the earliest attempts 
in this direction is connected with the classical Hyers-Ulam stability the- 
orem for e-convex functions. The most recent results are collected in |12j . 
Exact calculations with epsilons and sharp estimates are sometimes bulky 
and slightly mysterious. Some alternatives are suggested by actual infini- 
ties, which is illustrated with the conception of infinitesimal optimality. 

Assume given a convex operator / : X — > E U +00 and a point x in 
the effective domain dom(/) := {x £ X : f(x) < +00} of /. Given e > 
in the positive cone E + of E, by the e- subdifferential of / at x we mean 
the set 

d £ f(x) := {T S L{X,E) : (Vx G X)(Tx - Fx < Tx- fx + e)}, 

with L(X, E) standing as usual for the space of linear operators from X 
to E. 

Distinguish some downward-filtered subset § of E that is composed 
of positive elements. Assuming E and <§ standard, define the monad 
of S as [i{t§) := n{[0i £ ] : £ *= 0( *}- The members of /J.(<^) are positive 
infinitesimals with respect to $ . As usual, °S denotes the external set 
of all standard members of E, the standard part of S '. 

We will agree that the monad /i(<£") is an external cone over °i and, 
moreover, [i(S) n °E = 0. In application, S is usually the filter of order- 
units of E. The relation of infinite proximity or infinite closeness between 
the members of E is introduced as follows: 

e\ Ps e 2 <-> ei — e 2 € n{£) A e 2 — ei € /x(<f ). 

Since 

f| 5 e /(a)= |J dj(x); 

therefore, the external set on both sides is the so-called infinitesimal subd- 
ifferential of / at x. We denote this set by Df(x). The elements of Df(x) 
are infinitesimal subgradients of / at x. If the zero oiperator is an infin- 
itesimal subgradicnt of / at x then x is called an infinitesimal minimum 
point of /. We abstain from indicating § explicitly since this leads to no 
confusion. 
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Theorem. Let fx : X xY -> EU +00 and / 2 :FxZ->£U +00 be 
convex operators. Suppose that the convolution f 2 A j\ is inGnitesimally 
exact at some point (x,y,z); i.e., (f 2 A fi)(x,y) » fi{x,y) + f 2 (y,z). If, 
moreover, the convex sets epi(/i, Z) and epi(X, f 2 ) are in general position 
then 

D(f 2 A f 1 )(x,y) = Df 2 (y,z)oDf 1 (x,y). 
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Interval analysis for algorithms of idempotent and 
tropical mathematics^ 

Grigory L. Litvinov 

The idempotent interval analysis appears to be best suited for treating 
problems with order-preserving transformations of input data [1, 2]. It 
gives exact interval solutions to optimization problems with interval un- 
certainties in input data without any conditions of smallness on uncer- 
tainty intervals. Our aim to generalize results presented in [1, 2] for a very 

1 The work has been supported by the joint RFBR/CNRS grant 05-01-02807 and 
by the RFBR grant 05-01-00824. 



<)<S 



Grigory L. Litvinov 



general case of arbitrary algoritms of idempotent mathematics (in partic- 
ular, tropical mathematics) and algorithms over positive semirings (the 
scmifield of all nonnegative real numbers with usual operations is a typi- 
cal positive semiring) . Algorithms of this type are generated by a collec- 
tion of basic semiring/semifield operations, the well known star-operations 
x x* = I (& x (B x 2 ® x 3 (& . . . , and trivial operations. 

Theorem. Every algorithm of idempotent mathematics (and every 
algorithm over positive semirings) has an interval version. The complexity 
of this interval version coincides with the complexity of the initial algo- 
rithm. The interval version of the algorithm gives exact interval estimates 
for the corresponding output data. 

See [1, 2] for some examples. 

Note that for the traditional interval analysis the situation is oppo- 
site. For example, basic algorithms of the traditional linear algebra are 
plynomial but the corresponding interval versions are NP-hard and inter- 
val estimates are not exact. 
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1. The Maslov dequantization 

Let R and C be the fields of real and complex numbers. The well-known 
max-plus algebra R ma x = RU {—00} is defined by the operations x © y = 
max{x, y} and x y = x + y. 

The max-plus algebra can be treated as a result of the Maslov dequan- 
tization of the semifield R+ of all nonnegative numbers, see, e.g., [1,2]. 
The change of variables 

(1.1) x 1 — ► u — h log x , 

where h > 0, defines a map $h ■ R+ — > RU {—00}. Let the addition and 
multiplication operations be mapped from R + to RU {—00} by <f>;j, i.e. 
let 

u (Bh v — h\og(exp(u/h) + exp(v/h)), u v = u + v, 

= -00 = $ fc (0), 1 = = ^(1). 

It can easily be checked that u (Bh v — > max{ii, v} as h — > 0. Thus we get 
the semifield M max (i-e. the max-plus algebra) with zero = — 00 and unit 
1 = as a result of this deformation of the algebraic structure in K + . 

The semifield M ma x is a typical example of an idempotent semiring; 
this is a semiring with idempotent addition, 1.6. j x £B x — x for arbitrary 
element x of this semiring, see, e.g., [3-5]. 

The analogy with quantization is obvious; the parameter h plays the 
role of the Planck constant [2]. The map x 1— ► |x| and the Maslov dequanti- 
zation for H._|_ give us a natural passage from the field C (or M) to the max- 
plus algebra ]R max - Following [4], we will also call this passage the Maslov 
dequantization. In fact the Maslov dequantization is the usual Schrodinger 
dequantization but for imaginary values of the Planck constant (see, e.g., 
[4]). The passage from numerical fields to the max-plus algebra K max (or 
similar scmificlds) in mathematical constructions and results generates the 
so called tropical mathematics. The so-called idempotent dequantization is 
a generalization of the Maslov dequantization; idempotent dequantization 
generates the so-called idempotent mathematics, see, e.g. [4] for details. 



lr This work has been supported by the RFBR grant 05-01-00824 and the joint 
RFBR/CNRS grant 05-01-02807. 
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2. The dequantization transform 

This transform is denned in [6] . 

Let X be a topological space. For functions f(x) denned on X we 
shall say that a certain property is valid almost everywhere (a.e.) if it 
is valid for all elements x of an open dense subset of X. Suppose X is 
C" or R"; denote by R\ the set x = {(x u ...,x n ) G X | > for 
% = 1,2, ... ,n. For x = {x\, . . . , x n ) G X we set exp(x) = (exp(xi), . . . , 
exp(x n )); so if x G R™, then exp(x) G R". 

Denote by F(C n ) the set of all functions defined and continuous on 
an open dense subset U C C™ such that U D R™ . It is clear that ^C) 
is a ring (and an algebra over C) with respect to the usual addition and 
multiplications of functions. 

For / G J-(C n ) let us define the function fa by the following formula: 

(2.1) AOr)=/ilog|/(exp(x//i))|, 

where h is a (small) real positive parameter and x € R n . Set 

(2.2) fax) = lim fa(x), 

h— H-0 

if the right-hand part of ( |2.2| exists almost everywhere. We shall say that 
the function f(x) is a dequantization of the function /(x) and the map 
f(x) i — ► /(a:) is a dequantization transform. By construction, fa(x) and 
/(a;) can be treated as functions taking their values in K max - Note that 
in fact fa(x) and fax) depend on the restriction of / to K™ only; so in 
fact the dequantization transform is constructed for functions defined on 
M" only. It is clear that the dequantization transform is generated by the 
Maslov dequantization and the map i>-> \x\. Of course, similar definitions 
can be given for functions defined on K™ and Wl_. 

Denote by df the subdifferential of the function / at the origin. 

It is well known that all the convex compact subsets in R" form an 
idempotent semiring S with respect to the Minkowski operations: for 
A, B € S the sum A B is the convex hull of the union A U B; the 
product A B is defined in the following way: A&B — {x\x = a + b, 
where a <E A, b G B. In fact S is an idempotent linear space over R max 
(see, e.g., [4]). Of course, the Newton polytopes in V form a subsemiring 
TV in S. If /, g are polynomials, then d(fg) — df dg; moreover, if / and 

g are "in general position", then d(f + g) = 9/0 dg. For the semiring 
of all polynomials with nonnegative coefficients the dequantization trans- 
form is a homomorphism of this "traditional" semiring to the idempotent 
semiring N ' . 
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Theorem 2.1. If f is a polynomial, then the subdifferential df of f 
at the origin coincides with the Newton polytope of f. For the semiring 
of polynomials with nonnegative coefficients, the transform f i— ► df is a 
homomorphism of this semiring to the semiring of convex polytopes with 
respect to the well-known Minkowski operations. 

Using the dequantization transform it is possible to generalize this re- 
sult to a wide class of functions and convex sets, see [6] . Another approach 
based on complex analysis is due to A. Rashkovskii, see, e.g., [7,8]. 

3. Dequantization of linear operators and semigroups of linear 
operators 

The dcquantization transform can be rewritten in the following form: 
/ f(x) = lim Mog(| f(exp(x/h) |) = lim (1/s) • log(| /(exp(sx) |), 

where x G K™, and h, s — 1/h are real positive parameters. 

Our aim is to apply the dequantization transform to matrix elements 
of operator semigroups generated by linear operators. 

Suppose that S is a semigroup and s h- > n s is a linear representation 
of S in a complete (or quasicomplete) barreled locally convex space (by 
continuous operators). Denote by V the dual space to V and by (v',v) 
the value of a functional v G V on an element v G V. If s i — ^ n v >, v (s) = 
(v',tt s v) is a matrix element of n, then its dequantization 9 V >. V is defined 
by the formula: 

TT V ,. V = lim s ^oo(l/s) -log(| (v',TT s v) |). 

We discuss the cases S = M.+ or S = Z + . If v and v' are fixed, then 
tt v ', v eSU {oo}. 

Proposition 3.1. Let A be a linear operator in V, tt s = exp(sA), and 
dimV < oo. Then the set of all dequantizations {tt v ' iV } coincides with the 
set of real parts of all eigenvalues of A. 

There are generalizations of this result for the case dimF = oo. 

Suppose that for every v' G V there exists a number r > such that 
the set {r _s • (| (v',tt s v) \),s G S} is bounded for every v G V and the 
number r does not depend on v' G V. Then the representation it is called 
exponential. Note, that if V is a Banach space and ir is weakly continuous, 
then n is exponential. In the general case the spectral radius p„ of ir is 
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defined by the formula: 

Ptj = inf{r | r~ s n s v — > weakly for every v <E V as s — > +00} 

Proposition 3.2. 7/^4 is a bounded linear operator in a Banach space 
V , S = Z + , 7T = A s , then p„ = p(A) = lim^oo ||A S || 1 / S , i.e. p n is the 
traditional spectral radius of A. 

Theorem 3.1. If -k is exponential, then 

\ogp v = sup{w v , iV I v' e V',v e V}. 

Theorem 3.2. Suppose that A is a compact operator and ir s = A s , 
where s e 5 = Z + . T/ien £/ie se£ {??„',„} 0/ all dequantizations of ir 
coincides with the set of all numbers of the form log(| A |), where A runs 
the spectrum of A. 

4. Dequantization of set functions on metric spaces 

Let M be a metric space, S its arbitrary subset with a compact closure. 
It is well-known that a Euclidean d-dimensional ball B p of radius p has 
volume 

vo1 ^) ^ fTTrS) ^ 

where dis a natural parameter. By means of this formula it is possible to 
define a volume of B p for any real d [9]. Cover S* by a finite number of 
balls of radii p m . Set 

v d (S) := lim inf Vvol d (B Pm ). 

P^O p m <P Z ' 

m 

Then there exists a number Z) such that Ud(5) = for d > D and Vd{S) = 
00 for d < D. This number D is called the Hausdorff-Besicovich dimension 
(or HB-dimension) of S 1 [9] . Note that a set of non-integral HB-dimension 
is called a fractal in the sense of B. Mandelbrot. 

Denote by Af p (S) the minimal number of balls of radius p covering S. 
Then 

where D(S) is the HB-dimension of S. Set p = e~ s , then 
D(S) = hm^+coa/s) • log^ exp( _ s) (S). 

So the HB-dimension D(S) can be treated as a result of a dequantization 
of the set function Af p (S). 
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Let fi be a set function on M (e.g., a probability measure) and suppose 
that n{B p ) < oo for every ball B p . Let B xp be a ball of radius p having 
the point x <G M as its center. Then define fi x (p) := fj,(B XiP ) and let 

:= lim_ +00 - (1/s) • log(|^( e - s )|). 

This number could be treated as a dimension of M at the point x with 
respect to the set function jj,. There are many dequantization procedures 
of this type in different mathematical areas. In particular, V.P. Maslov's 
negative dimension [10] can be treated similarly. 

5. Dequantization of the Fourier-Laplace transform 

It was noticed by V.P. Maslov (see, e.g., [1-4]) that the Lcgcndre (or 
Legendrc-Fenchel) transform can be treated as an idempotcnt (or trop- 
ical) version of the Fourier-Laplace transform. It seems to be interesting 
to note that the Legendre transform can be constructed from the Fourier- 
Legendre transform directly by means of the Maslov dequantization. 

6. Dequantization of geometry 

An idempotcnt version of real algebraic geometry was discovered in the 
report of O. Viro for the Barcelona Congress [11]. Starting from the idem- 
potent correspondence principle [2], O. Viro constructed a piecewise-linear 
geometry of polyhedra of a special kind in finite dimensional Euclidean 
spaces as a result of the Maslov dequantization of real algebraic geometry. 
He indicated important applications in real algebraic geometry (e.g., in the 
framework of Hilbcrt's 16th problems) and relations to complex algebraic 
geometry and amoebas in the sense of I. M. Gelfand, M. M. Kapranov, and 
A. V. Zelevinsky. Then complex algebraic geometry was dequantized by 
G. Mikhalkin and the result turned out to be the same; now the new geom- 
etry is called tropical algebraic geometry. In particular, tropical varieties 
are results of a dequantization procedure (generated by the Maslov de- 
quantization) applied to algebraic varieties. There are many applications, 
see, e.g., [11-13,5]. 

7. Remark 

It would be nice to find new dequantization procedures related to the 
Maslov dequantization. 
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